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We argue for the existence of plasma-baUs - meta-stable, nearly homogeneous lumps of 
gluon plasma at just above the deconfinement energy density - in a class of large 
^ ' confining gauge theories that undergo first order deconfinement transitions. Plasma-balls 
decay over a time scale of order A^ by thermally radiating hadrons at the deconfinement 
temperature. In gauge theories that have a dual description that is well approximated by 
a theory of gravity in a warped geometry, we propose that plasma-balls map to a family 
of classically stable finite energy black holes localized in the IR. We present a conjecture 
^ ■ for the qualitative nature of large mass black holes in such backgrounds, and numerically 
construct these black holes in a particular class of warped geometries. These black holes 
have novel properties; in particular their temperature approaches a nonzero constant value 
at large mass. Black holes dual to plasma-balls shrink as they decay by Hawking radiation; 
towards the end of this process they resemble ten dimensional Schwarzschild black holes, 
which we propose are dual to small plasma-balls. Our work may find practical applications 
in the study of the physics of localized black holes from a dual viewpoint. 
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1. Introduction and Summary 

In this paper we study the spectrum of localized excitations in confining large N 
gauge theories with energies of order A^^, the spectrum of static black holes in warped 
geometries with an effective infrared wall, and the connection between these two spectra, 
using generalizations of the AdS/CFT correspondence. In the rest of this section we will 
first introduce each of these questions and then explain the relationship between them. 
Readers who are only interested in gauge theory aspects may read sections 1.1, 2 and 3 
and appendices A, B, C, which are independent of the rest of the paper. 

1.1. Meta-stable Bubbles of Plasma in Confining Large N Gauge Theories 

It has long been believed [|ll that confining large A^ SU (N) gauge theories are dual to 
weakly coupled string theories with a string coupling Qs oc 1/A^. The string duals possess an 
exponential tower of long lived excitations (perturbative string states) which are identified 
with glueballs. However, at least some perturbative string theories also possess other long 
lived excitations - black holes - at energies of order l/g"^- In this paper we argue that such 
configurations have analogues in a class of confining large A^ gauge theories (for earlier 
related works see 

Consider a confining large A^ gauge theory, with mass gap Agap, whose thermal decon- 
finement phase transition is of first order (e.g. pure Yang Mills theory at large A^ 0,3)- In 
this paper we argue that, provided the theory in question obeys one additional condition 
(see section 2), it hosts long lived excitations at all masses m ^ N'^Agap. We call these 
meta-stable configurations plasma-balls. Plasma-balls are spherical, approximately homo- 
geneous lumps of deconfined plasma fiuid, whose energy density is just above the critical 
density (which is the energy density at the deconfinement phase transition temperature 
Td] this energy density is of order A^^ in the large A^ limit). These lumps are static because 
the pressure of the plasma in such theories vanishes at the critical energy density. 

Of course a plasma-ball is not completely stable; eventually it decays into a gas of 
hadrons. In section 3 we study the A^ dependence of the processes that contribute to this 
decay, and conclude that the hadronization of the plasma-ball is very slow in the large 
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limit. Even though the density of the plasma-ball is of order iV^, and the number of 
gluon-gluon interactions per unit time within the plasma-ball is of order A^^, glueballs are 
produced from these collisions at a rate that is independent of A^. Intuitively, only 
of the gluon-gluon collisions - collisions between a gluon and its color anti-partner - can 
form a singlet state that can escape from the plasma. 

Large counting is insufficient to determine the precise rate of radiation of glueballs 
as a function of (for instance) their mass; however the thermal nature of the plasma-ball 
suggests that these rates are controlled by the Boltzmann factor at the temperature of the 
plasma-ball. We conjecture that this is indeed the case. In a class of examples described 
below, the dual gravitational description of plasma-balls confirms this conjecture. 

We have already explained that the rate of radiation of any given glueball from the sur- 
face of a plasma-ball is independent of A^. Using the conjecture of the previous paragraph 
and the expectation that the density of glueball species at high energies is proportional to 
exp(£'/Tff) where Th is the Hagedorn temperature, the rate of loss of energy integrated 
over all glueball species is proportional to / exp {jy^ ~ 't^) ^-^^ deconfinement 
temperature is strictly smaller than the Hagedorn temperature in theories that undergo 
first order deconfinement transitions [^,0, this is a finite number of order A^*^. As a conse- 
quence, the lifetime of a plasma-ball of radius R {R^ 1/ Agap) is of order N^R. 

Note that QCD at zero chemical potential (for baryon number) does not undergo a 
first order phase transition as a function of temperature. Thus, we do not expect plasma- 
ball-like meta-stable configurations to be created in experiments like RHld. See appendix 
B.5 for further discussion. 

1.2. Black Holes in the IE 

We start this subsection by characterizing the class of warped geometries whose black 
holes we study. Consider any solution to Einstein's equations coupled to appropriate 
matter fields, whose metric can be put in the form 

ds^ = a'L^ {W'^{u)dxl + ds^) , (1-1) 

^ Dynamical plasma configurations which are dual to dynamical black holes could perhaps be 
related to RHIC, as suggested by Nastase and collaborators [^. In this context note that the decay 
products of the RHIC fireball have a distribution which is close to thermal, at approximately the 
deconfinement crossover temperature. 
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where fi = 0,1,2, ■■■ ,p, a'L"^ is a constant with units of length squared, and (isf^^ is the 
metric on an internal manifold, one of whose coordinates is u, and whose constant u slices 
are compact. The variable u in ( |1.1D has the range uq < u < oo and the function W{u) 
increases monotonically from a positive nonzero value at uq to infinity at w = oo. We 
assume that ( p. . 1|) is smooth and without boundaries everywhere, including at the IR wall 
■u = uo & We further assume that the space in question has well defined thermodynamics 
(in particular it admits the definition of an ADM mass), and that the spectrum of fiuctu- 
ations about (|1.1[) is gapped. Examples include the confining backgrounds constructed in 
110 • i 



We are interested in backgrounds ( p. . 1|) which have an additional property; they host 
black brane solutions of finite energy density @ at all energy densities larger than a critical 
value pe- Moreover, these branes are required to be stable (and have positive specific heat) 
above ps > Pe, and to have negative free energy density above Pc > Ps- Let Td denote the 
temperature of the black brane of energy density pc- Assuming that there are no other 
phases of the theory, it follows that the thermodynamics in this background is dominated 
by a graviton gas about ( |1.1D for T < Td, and by the black brane for T > Td; the system 
undergoes a first order phase transition at T = T^. 

In this paper we conjecture that spaces that obey all these properties possess a one 
parameter set of spherically symmetric (in p dimensions) black hole solution^ labeled by 
their mass, and that, in the large mass limit, these black holes have the following properties: 

1. Their volume in p dimensions is proportional to their mass, i.e. the radius of the black 
hole in p dimensions scales with mass like {m/ pc)^^^. 

2. In the interior (meaning away from the edges in the IR^ spatial directions, not the 
black hole interior behind the horizon) these black hole solutions approximate the 
black brane at energy density pc- 



This is possible because a fc-cycle of the internal manifold shrinks to zero size at u = ■uo, so 
that, locally, u — uq may simply be thought of as the radial coordinate of an IR*^^^ component of 
the geometry. 

^ The background of actually does not have a mass gap |11], but we expect that our analysis 
should still apply to it. 

^ More precisely, there exist smooth gravitational solutions with horizons, asymptoting to ( p..l[ ), 
which preserve translational invariance in IR^ and have a finite energy density above that of the 
background (|1 . 1|) . 

^ Black holes in this context were first discussed in Q. 
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3. In the vicinity of their edge the solutions reduce to a domain waU in the directions, 
that interpolates between the black brane at the critical energy density pc at one end, 
and the background (|n|) at the other end. 

We emphasize that the black holes described in the previous paragraph have ther- 
modynamical properties that are qualitatively different from black holes in flat space; in 
particular their temperature tends to a finite value in the limit of infinite mass. 

In section 4 we use intuitive arguments to motivate the conjecture described above. 
In addition, as concrete evidence for our conjecture, in section 8 we present a numerical 
construction of the domain wall about a particularly simple background of the form ( |1 . 1|) . 
The background @ we study is a solution to Einstein's equations with a negative cos- 
mological constant but no other matter fields, and asymptotes to a Scherk-Schwarz circle 
compactification of Lorentzian AdSd+2, or (at finite temperature in Euclidean space) to 
a torus compactification of Euclidean AdSd+2- An infinite series of translationally invari- 
ant (in d — 1 dimensions) solutions with these asymptotics is known. The topology of 
each of these solutions is IR'^ times the solid torus; the distinct solutions are labeled by 
which cycle of the boundary torus is 'filled in' in the solution. When the boundary torus 
is rectangular, the thermodynamics is dominated by the solution that fills the smaller 
of the (1,0) and (0,1) cycles. A phase transition between these two solutions occurs 
when the boundary torus is a square. In section 8 we give an overview of the numerical 
construction of the domain wall solution that interpolates between these two solutions 
at the phase transition temperature, and compute the (positive) surface tension of this 
domain wall. More details of our construction may be found in appendices F and G, 
and an example of the program used to generate the solutions may be downloaded from 



http : / / schwinger . harvard . edu/ ^wiseman/IRblackholes/ . 



1.3. Plasma-Balls as Black Holes 

Dual string descriptions of large SU{N) gauge theories are expected to exist on 
general grounds from 't Hooft's analysis of the large limit [Q. The dual string coupling 
is always proportional to 1/A^, so the interactions between strings are small in the large 
limit. Unfortunately the string duals to most familiar confining gauge theories (such 
as pure SU{N) Yang-Mills theory) have yet to be discovered; furthermore, the relevant 
dual dynamics may be argued to be controlled by a strongly coupled worldsheet theory 
that does not admit a truncation to the lowest (gravitational) sector. In several examples 
H,p!2|,|9| JT0|] , however, it turns out to be possible to modify standard confining gauge theories 
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by coupling them to a judiciously chosen set of additional degrees of freedom. The energy 
scale of 'new physics' is a free parameter, which is most conveniently labeled by the value 
of the 't Hooft coupling A = Qym-^ scale of new physics. In the limit A ^ the 

scale of new physics is much higher than the mass gap, and, at low energies compared to 
the scale of the new physics, we recover the (strongly coupled) dynamics of the original 
confining theory. In the opposite limit of large A, the new degrees of freedom give large 
anomalous dimensions to all operators that create glueballs dual to string oscillators. As a 
consequence such glueballs have masses that are parametrically higher than the graviton, 
and the low energy dual dynamics reduces to a theory of gravity (usually ten dimensional 
supergravity) in a warped background. In specific examples |]8|,p!2|,|9|,p!0[] this gravitational 
dual may explicitly be identified using a generalization of the AdS/CFT correspondence. 
In each case the resulting gravitational theory lives in a warped geometry of the form 
( p. . 1|) (for a p + 1-dimensional confining gauge theory) that obeys all the conditions of the 
previous subsection, and so, according to our conjecture, hosts localized black holes. Since 
the 'black brane nucleation' phase transitions in these geometries may be identified with 
deconfinement transitions in the dual field theory [0], it follows that these localized black 
holes are dual to plasma-balls (see P,p!B|,P| for related remarks). 




Evapouraies lo glueballs at surface 





Figure 1: The plasma-ball and its decay via hadronization (left), and the 
localized black hole and its decay via Hawking radiation (right). 
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The decay of the plasma-baU by hadronization maps to the decay of its dual black 
hole by Hawking radiation. It follows that, at least at large A, the decay of plasma-balls 
is a thermal process at the temperature T^, in agreement with the guess described in §1.1. 

Even in warped backgrounds, high energy graviton-graviton collisions are likely to be 
dominated by black hole formation It follows that, in confining theories at large A, 

high energy hadron-hadron collisions are dominated by plasma-ball production. Of course, 
plasma-ball formation does not dominate high energy glueball-glueball collisions at small 
A. In Section 6 we trace these different behaviors to qualitative differences in the hadron 
parton distribution functions at small and large A [|n[ . 

The duality with plasma-balls could yield interesting lessons for black hole physics. As 
we have explained above, at large A the production of a plasma-ball in high energy hadron- 
hadron collisions, and its subsequent decay by hadronization, is dual to the production of 
a black hole in high energy graviton-graviton scattering, and its subsequent decay by 
Hawking radiation. As the first of these processes is manifestly unitary, its dual must be 
as well. It is possible that this connection could be pursued further to draw more detailed 
lessons about the nature of the Hawking radiation process. 

As a second application, we argue in section 7 that one of the most striking features 
of black hole physics - the universally absorptive nature of black holes - is crucially tied 
to a feature of the hadronic parton distribution function at large A that is absent at small 
A. It follows that while black hole like configurations continue to exist in at least some 
warped geometries with curvatures of order the string scale (those that are dual to gauge 
theories at small A) , they have qualitatively different dynamics from their large A cousins. 
In particular they are no longer black. 



2. The Plasma-ball as a Stable Lump of Plasma Fluid 

In this elementary sub-section we will argue for the existence of a meta-stable localized 
lump of plasma fluid (in the deconflned phase) - a plasma-ball - in certain large theories 
that undergo flrst order phase transitions. More details may be found in appendix A. 

Consider an isolated spherical ball of static plasma fluid of radius R ^ 1/Agap in 
p spatial dimensions. In equilibrium its effective temperature T and the pressure P are 

7 



both uniform in the bulk of the bali. In order for the ball to be static, the pressure of 
the plasma in the ball must precisely balance the tension of the domain wall separating 
the two phases (the "surface tension" of the plasma fluid; see appendix B.l for a precise 
deflnition). Denoting this tension by E, the forces balance when 

u:^_^W-^T. = ^^RP-^P ^ P= (2.1) 
p — 1 R 

where ujp-2 is the surface area of the unit {p — 2)-sphere, and ujp-2/{p — 1) is the volume 
of the unit {p — l)-ball. Since the bubble is in the deconfined phase, its pressure and 
surface tension are both of order A^^ 0, and are functions of the plasma-ball temperature. 
According to ( p.l| ), static plasma-balls exist at asymptotically large R if and only if the 
pressure of the deconflned phase vanishes at a finite temperature. 

While the precise functional form of P{T) (the pressure of the deconfined phase as a 
function of temperature) depends on details, this function is constrained by thermodynamic 
considerations. First, a simple thermo dynamical identity (see appendix A) ensures that 
the pressure increases monotonically with temperature, provided that the specific heat 
of the deconfined phase is positive. Second, as pressure is continuous across a phase 
transition, the pressure of the deconfined phase is 0{1) (rather than the generic 0{N'^)) at 
the deconfinement temperature T^. Provided that the phase transition is of first order (an 
assumption we will make in most of the rest of this paper) the specific heat, and so dP/dT, 
are positive and of order A^^ at T = (see appendix A) 0. It follows that the pressure 
of the deconfined plasma vanishes at T = T^^ — 0{1/N'^). We conclude that uniform, 
asymptotically large lumps of fiuid plasma are static at the deconfinement temperature 
in large gauge theories that undergo first order phase transitions. Provided that E 
is positive, large but finite static spherical lumps of plasma fiuid exist at a temperature 

^ In this section we ignore energy loss by 'radiation' from the surface of the plasma-ball. In the 
next section we will argue that this approximation is justified in the large N limit. The energy 
loss by radiation to the bulk is 0(1); as the energy density within the plasma ball is 0{N'^), the 
temperature gradients induced by surface radiation are negligible in the large N limit. 

As we review in appendix A, the pressure is simply minus the free energy density of the 
bubble. The fact that the surface tension should scale as N'^ is predicted from general arguments, 
since planar diagrams should contribute to it, but it has not yet been directly verified by lattice 
computations [Q. 

More precisely this is true of the limit of the specific heat, as T approaches Td from above. 
This limit determines the speed of sound of the deconfined plasma fluid within the plasma-ball. 
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slightly above the deconfinement temperature, and have an energy density slightly above 
the critical energy density. 

These static lumps are also hydrodynamically stable 0, provided that the effective 
surface tension is positive. As we have assumed the deconfinement transition to be of first 
order, the deconfined phase continues to exist (as a meta-stable phase) at temperatures 
below Td- As we have explained in the previous paragraph, the pressure of the deconfined 
phase is positive for T > and negative for T < T^; it follows that plasma-balls are 
stable against homogeneous expansion or contraction. Local stability of the plasma fiuid 
against density fiuctuations in the bulk of the ball follows from the positivity of the speed 
of sound at the phase transition temperature (see appendix A). Linearized long wavelength 
fiuctuations of the surface also obey a wave equation; the squared speed of sound of these 
surface waves is proportional to the effective surface tension; as a consequence the plasma- 
ball is stable to fiuctuations of the boundary provided that this surface tension is positive 
and approaches a finite value at the deconfinement temperature. 

Note that according to (|2.1J ) the pressure (hence temperature) of a plasma-ball with 
positive surface tension is a decreasing function of its radius (hence mass) , so plasma-balls 
have negative specific heat 0. Of course, plasma-balls are not completely stable; they 
eventually hadronize. In the next section we discuss this hadronization process. Additional 
discussions of various aspects of plasma dynamics may be found in appendix B. 

In this section we have, so far, discussed plasma-balls that carry no conserved charges 
besides their energy, and we will focus on such configurations in the remainder of the paper. 
However the hydrodynamical construction of stationary lumps of plasma generalizes in 
a straightforward manner upon adding other conserved charges. As a simple example 
consider a radially symmetric lump of plasma rotating with angular velocity to (about the 
origin) in p = 2 spatial dimensions. Force balance yields the equation 

^ = p{r)uj\ (2.2) 

They resemble the bubbles of the deconfined phase which appear when the energy density 
of the system is gradually raised through the first order phase transition, except for the fact that 
they live in the vacuum. 

Note that the (negative) specific heat of the plasma-ball as an object is distinct from the (pos- 
itive) specific heat per unit volume of the deconfined phase of which the plasma-ball is composed. 
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where P{r) and p(r) are the pressure and density at radial distance r, and we have used 
the non-relativistic approximation appropriate for smaU rotational velocities. Substituting 
p as a function of P using the equation of state, we find 



This may be regarded as an equation for P{r), and potentially it can have two qualitatively 
distinct classes of solutions. In solutions of the first class the variable r has the range 
(0, R). Such a solution may be thought of as a rotating plasma-ball. In the second class 
of solutions r has the range and ( |2.3| ) is to be solved subject to the boundary 

condition P{Ri) = —T,/Ri (we assume i?, Ri ^ 1/ Agap). Such a solution is best thought 
of as a rotating plasma-ring, filling an annulus in the plane. 

3. Hadronization of Plasma-Balls at Large 

In this section we will compute the N dependence of the rate at which a plasma- 
ball hadronizes (decays into glueballs). We model the plasma-ball as a gas of initially 
thermally distributed uncorrelated gluons, whose collisions sometimes produce glueballs. 
We then determine the dependence of all Feynman graphs (with arbitrary numbers of 
interaction vertices) that contribute to glueball production. We find that the rate at which 
any particular glueball is produced is independent of N. 

Of course, the gluonic constituents of a plasma-ball are not really uncorrelated; how- 
ever all correlations may be switched off without encountering a phase transition (for 
instance by going to high temperatures in an asymptotically free theory), so they should 
not affect the A^-scaling of the gluon production rates computed in this section. The 
arguments of this subsection are valid at all orders in perturbation theory; however, ex- 
perience with counting powers of in the 't Hooft limit (e.g. in exactly solvable matrix 
models) suggests that the powers of obtained from this argument will be correct even 
non-perturbatively. In specific examples we will provide evidence that this is indeed the 
case in sections 5-8. 
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3.1. Counting Powers of N 

Consider a connected Feynman graph that describes n initial gluons scattering into 
m final gluons and k glueballs. The graph in question has n + m external gluon lines, and 
includes k insertions of glueball creation operators such as TrF^^F^^^. These operators are 
normalized so that their two point functions are of order A^'^ to ensure that they create 
glueballs with unit probability; this means that each such operator insertion in a Feynman 
diagram appears with an extra factor of 1/N compared to an insertion of an interaction 
vertex. The contribution of this graph to the inclusive probability for glueball production 
is obtained by squaring the graph and summing over all initial and final gluon states. We 
will now determine the N dependence of the result of this process. 




Figure 2: A typical sewn graph contributing to the glueball production 
rate, in double-line notation. 

Consider the graph in question, drawn using the 't Hooft double line notation, together 
with its CPT conjugate (a graph in which fundamental and anti-fundamental indices are 
interchanged). Sew these two graphs together, as in figure 2, by attaching every free gluon 
line in the first graph to the corresponding conjugate gluon line in the conjugate graph; 
this gluing preserves the flow of all color indices. The resulting graph has no free gluon 
lines, so the faces may be filled in to form a genus g Riemann surface in the usual manner. 
The N dependence of the contribution of this graph to inclusive glueball production is 
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obtained by freely summing over all indices of the sewn graphed, and so is proportional to 
jY2-2g-2fc (^^gjj^g standard 't Hooft counting and the normalization of the glueball creation 
operators). The leading behavior at large is given by planar graphs (^r = 0) and is 
proportional to 

iV2-2^ (3.1) 

3.2. k=0: Gluon Mean Free Time 

According to ( |3.1D the effective 'number' of gluon-gluon collisions per unit time, in a 
large plasma-ball of volume V , is of order A^^. It follows that the rate at which any given 
gluon undergoes collisions is of order A^'^ (recall that the plasma-ball contains of order A^^ 
gluons). As a consequence, the relaxation time scale of the plasma fluid is of order A^*^. In 
appendix C we verify this result at lowest order in perturbation theory. 

3.3. k=l: Glueball Production Rate 

The plasma-ball radiates glueballs by one of two mechanisms. Every once in a while 
a gluon with indices {i,i) collides with a gluon with indices (j, i) near the surface of the 
plasma-ball; this collision can produce a glueball, which then escapes from the plasma-ball. 
Alternatively, occasionally an energetic gluon shoots out of the surface of the plasma-ball; 
the string that attaches it to the plasma-ball can then snap, allowing it to escape into 
the bulk as a glueball. Glueball production from each of these processes is governed 
by the Feynman graphs studied in §3.1 above. According to (|3.1| ), in the large A^ limit 
glueball production is dominated by graphs with /c = 1, and takes place at a rate that 
is independent of A^. In appendix C we study a simple Feynman graph that contributes 
to glueball production, to verify this conclusion at weak coupling (the general analysis is 
valid at any value of the 't Hooft coupling). 

In this subsection we have focused on the radiation of glueballs from plasma-balls. 
However the arguments of this section apply equally well to the time reversed process. 
Consider a glueball incident upon the plasma ball from outside. Equation ( pT|) implies 

This may be argued as follows. Color indices attached to the m final gluons must be summed 
over in summing over all final states. Color indices attached to the n initial gluons must be summed 
over in summing over all initial gluons - in this step we use the fact that the plasma-ball has a 
finite density of gluons of any given color variety. All internal color indices must be summed over 
in squaring the original Feynman graph. 
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that the interaction cross section for this gluebaU, per unit length traversed through the 
plasma ball, is of order N^. We expect this interaction process to lead to the glueballs 
dissolving into the plasma-ball. Similarly, glueballs formed in a gluon collision far from the 
surface of a large plasma-ball will dissolve before they escape, so the hadronization rate is 
proportional to the surface area of the plasma-ball rather than its volume. 

3.4. k=2: Friction 

As we have discussed in the previous subsection, the cross section for a glueball travers- 
ing the plasma-ball to dissolve into the plasma is of unit order. Glueballs also interact with 
the plasma in a more elementary manner; gluons incident on the glueball knock it around, 
slowing down the glueball and causing it to jiggle around as it traverses the plasma-ball. 
Such interactions are governed by diagrams with one incoming and one outgoing glueball, 
have k = 2 and so are of order l/N"^. 



4. Localized Black Holes in Warped Backgrounds 

We now turn to a study of uncharged, static, finite energy (rather than finite energy 
density) black holes in the backgrounds ( p. . 1|) , that are localized in the IR region of the 
radial direction as well as in the IR^ directions, in the range of parameters where gravity in 
the background ( p..l|) is a good approximation to the theory. At low enough energies such 
solutions are well approximated by ten dimensional Schwarzschild black holes of radius 



Rs ^ V a'L, centered aX u = uq in (|1.1|) . What is the qualitative nature of the static black 



hole solutions at higher energies? A consideration of the opposite, high energy limit gives 
us a clue. The confining backgrounds ( |1 . Ij ) host infinite energy translationally invariant 
(in the field theory IR^ directions) black brane solutions, which undergo a first order phase 
transition at some finite temperature^. The thermodynamic properties of these black 
branes are clearly extensive (for instance, their energy is proportional to their volume 
in IR^). Recall also that the graviton (as well as all other fields) is massive in ( p. . 1|) ; as 
a consequence all correlation functions decay exponentially in this background. Putting 
these facts together, it is natural to expect that a black hole whose energy is very large, 
rather than strictly infinite, is a large spherical lump in IR^ which in its interior (meaning 



Note that the deconfinement transition is always of first order when the gravitational ap- 
proximation is valid, since it occurs far below the Hagedorn temperature which is related to the 
string tension. 
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away from its edges in the directions - not the black hole interior behind the horizon) 
closely resembles a translationally invariant black brane. Furthermore, taking the energy 
and the radius of the lump to infinity and focusing on the edge of the lump, we expect 
the solution to tend to a planar domain wall interpolating between the confining gravity 
solution (which is just (p. .11)) on one side and the deconfined homogeneous gravity solution 
(which is the black brane in ( p. . 1|) ) on the other side. 
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Figure 3: The localized black hole. The vertical axis is the radial coordi- 
nate, and the horizontal axis x is one of the spatial IR^ coordinates. 



Of course, uncharged homogeneous black brane solutions exist over a range of temper- 
atures. However, the force balance arguments of §2, which apply to warped backgrounds 
whenever they possess a well defined boundary stress tensor, tell us that such a large bub- 
ble can be static only if the relevant black brane has a pressure that vanishes in the large 
size limit. This is true for the black brane precisely at the deconfinement temperature (see 
§2 and appendix A). Indeed, the force balance equation ( |2.1|) is ensured by the conserva- 
tion of the (boundary) stress-energy tensor (see appendix B.l) which in turn is a direct 
result of diffeomorphism invariance of the bulk theory in warped gravitational backgrounds 
(see, for instance, [1T5| , p!^ and references therein, and see section 8 and appendices E and 
F for how this works in a specific example). As a consequence, in the classical gravity 
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approximation (i.e. to leading order in 1/A^), we expect to find large black holes whose 
interior locally reduces to a homogeneous black brane at the deconfinement temperature. 
The width of the "domain wall" separating the two regions is expected to be of the order 
of the curvature scale (which is also the scale of the mass gap in these theories). 

Recall that any system undergoing a first order phase transition goes through con- 
figurations that consist of bubbles of varying size of one phase inside the other. The 
pancake-like solutions we have described above are simply increasingly large bubbles of the 
deconfined phase within the confined phase, except that the confined phase is replaced by 
the vacuum. Within the classical supergravity approximation the thermal confined phase 
is indistinguishable from the vacuum (because its energy is 0{1) rather than 0{N^)), so 
this last replacement changes nothing. As in section 2, we expect these big black holes to 
be stable if their surface tension is positive. 

Motivated by these considerations, we conjecture that warped backgrounds of the form 
described in §1.2 always host large pancake-shaped black holes of the form described above 
(and in the introduction) and depicted in figure 3. In section 8 we will supply evidence for 
our conjecture by constructing the domain wall that interpolates between the vacuum and 
the black brane at the deconfinement temperature, in a particularly simple background. 
We also show that in this specific background the surface tension is positive, so we expect 
the corresponding black holes to be stable. 

Note that the black holes described in this subsection have a finite temperature 
in the large mass limit. Moreover, assuming that their surface tension is positive, they 
have negative specific heat for a rather prosaic reason (see ( |2.1J )); smaller black holes have 
to be hotter in order to balance the surface tension. Note that this behavior matches 
smoothly onto the negative specific heat of very small black holes (which should behave 
like Schwarzschild black holes, see §7.1). 

In this section we have focussed on a study of black holes with vanishing angular 
momentum (and other conserved charges). However, it seems likely that reasoning identical 
to that presented at the end of section 2 could be used to predict the detailed properties 
of rotating black holes and black rings (whose horizon contains an component, similar 



The domain walls constructed in section 8 do, indeed, have positive surface tension. We are 
not, however, aware of a rigorous argument that guarantees this will always be the case. If negative 
surface tension domain walls do exist, we would expect them to possess Gregory-Lafiamme like 
instabilities localized at the surface of the black hole, in analogy with the Gubser-Mitra conjecture 
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to the ones of |^), in the backgrounds studied in this section. We wiU leave a detailed 
analysis of such solutions to future work. Here we will merely note that the rotating 
plasma-ring like solutions described in section 2 exist only when p > 2 for the excellent 
reason that when p = 1 the plasma-fluid has no angular direction in which to rotate. As 
we have explained, gravitational backgrounds with p + 1 dimensional Lorentz invariant 
sections, of the form considered in this section, are at least p + 3 dimensional (the two 
additional dimensions refer to the radial u coordinate plus the compact part of the internal 
dimensions, which must be at least one dimensional since some dimensions must shrink in 
the IR) . It follows that an analysis of the sort presented in section 2 predicts black ring 
solutions only in 4 -|- 1 and higher dimensions. This agrees nicely with the fact that (at 
least in flat space) black rings do not exist in 3 -f- 1 dimensions, and it predicts that (at 
least in some gravitational theories) black rings should exist above 4 + 1 dimensions. 

Similarly, we expect to be able to generalize our solutions by adding charge (global 
charge in the gauge theory, which maps to a local charge in the gravitational dual) . Again 
we will leave such generalizations to future work. We also expect that there may be 



phenomenological applications for these localized solutions [21 



5. Localized Black Holes as Plasma-balls 

We now turn to a dual gravitational description of plasma-ball dynamics in bulk 
duals of conflning large gauge theories at large A. As mentioned in the introduction, 
there exist several examples of conflning gauge theories that have purely gravitational dual 
descriptions in the large A limit. Examples include (a speciflc UV completion of) 4 + 1 
dimensional super- Yang-Mills theory compactifled on a circle with anti-periodic boundary 
conditions for the fermions , "little string theories" (coming from type IIB NS5-branes) 
compactifled on an S"^ [|T^, and cascading gauge theories |j2^,0]. The geometries dual to 
these systems differ qualitatively in several respectslll; however they all share the general 
properties described in §1.2. According to the conjecture of the previous section, these 
backgrounds all possess localized black hole solutions, that closely resemble the deconflned 
black brane at their center. 

For example, in the dual of a 2 + 1 dimensional confining theory which we will discuss in 
detail in §8, the ring would be a solution of a 4 + 1 dimensional gravitational theory (with negative 
cosmological constant), with horizon topology x S'^. 

This is natural, as the UV dynamics of the various systems differs qualitatively. 
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A standard entry in the AdS/CFT dictionary identifies the translationaUy-invariant 
black brane at any given temperature with the deconfined phase of the dual gauge theory, 
at the same temperature [|]. It follows that the localized black holes referred to above are 
dual to localized lumps of plasma at approximately the deconfinement temperature, which 
are plasma-balls. 

In order to build a more detailed dictionary between bulk and boundary degrees of 
freedom, it is useful to contemplate a Gedanken experiment. Consider a black hole, of 
the form described in section 4, localized near the origin of p dimensional space. Imagine 
scattering a graviton in some mass eigenstate (about the original confining background) off 
the black hole at a given p dimensional momentum. The wave-function that describes this 
scattering process has four components; an incident wave, a reflected wave, a transmitted 
wave (the part of the wave that stays completely outside the horizon by going around the 
black hole, including a part going 'above' it in the radial direction), and a component of 
the wave absorbed into the black hole. Each of these components has a natural dual inter- 
pretation; a glueball incident on a plasma-ball may be reflected from the ball, transmitted 
through the ball or may be absorbed into the plasma-ball (thereby dissociating into glu- 
ons). We thus conclude that gravitons outside the horizon 'above' the black brane region 
of the localized black hole map to glueballs inside the plasma-ball 111 il . On the other 
hand, anything near the horizon (or any possible degrees of freedom behind the horizon) 
maps to gluon or deconfined degrees of freedom. 

As we have already argued in §3, lumps of the deconfined phase have zero overlap with 
glueballs only at infinite A^. At finite lumps of gluon plasma overlap with glueballs; 
this mixing explains the decay of plasma-balls by hadronization. These statements have a 
direct gravitational analogue; black holes at finite temperature mix with the gravitons at 
finite Qs, a phenomenon that explains the decay of non-extremal black holes by Hawking 
radiation. 

Quantum mechanically the momentum of such a graviton fluctuates. The lowest order dia- 
grams that contribute to such fluctuations are of order Qs, leading to a fluctuation probability of 
order ~ in agreement with the analysis of §3.4. 

Gravitons located far in the UV correspond to glueballs that are much smaller than A'g^p- 
Of course, such special glueball configurations have the propensity to expand, which is dual to 
the propensity of a graviton at large u to fall into the horizon. 
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6. Plasma-Ball Dynamics from Black Holes 

6.1. Plasma-Ball Hadronization as Hawking Radiation 

We have argued for the existence of stable pancake-lilce localized black hole solutions 
in backgrounds of the form ( p. . 1|) ; we have also argued that these black holes are dual to 
plasma-balls. Of course, these black holes are stable only classically; quantum mechanically 
they decay via Hawking radiation. 

Hawking radiation from the 'flat' surface of the pancake (see flgure 3), which goes 
out in the radial direction, is simply reflected back into the black hole (by the effective 
potential arising from the warp factor in the geometry of ( p. . 1|) , or, for special modes, by 
reflecting boundary conditions in the large u boundary of (|1 . 1| ) ) . This fact has a simple 
interpretation in the dual gauge theory; as discussed in §3, glueballs created by collisions 
between gluons in the bulk of the plasma-ball dissociate before they are able to escape out 
of the plasma-ball. The black hole loses energy only from Hawking radiation at its edge, 
reflecting the fact that only glueballs produced by gluon-gluon collisions near the surface 
of the plasma-ball can escape. 

From a gauge theory point of view, the decay of a plasma-ball is a non-perturbative 
process that is diflicult to study even numerically on the lattice. At large A the AdS/CFT 
correspondence maps this hadronization to the quantitatively well understood process of 
Hawking radiation, permitting a quantitative analysis of decay rates and branching ratios 
in any model of interest. In particular, this conflrms the thermal nature of the decay of 
the plasma-ball, at least at large A. 

In the fleld theory analysis adding a flnite number of quark flavors does not change 
the analysis, and the same is true also in the string theory duals. Adding flavors in the 
duals is done by adding D-branes, and in the 't Hooft large limit with a flnite number 
of D-branes, the D-branes have a very small back-reaction and do not affect any of the 
qualitative discussion above. Note that the Hawking radiation in this case is expected 
to be partly into closed string modes (glueballs), such as gravitons, and partly into open 
string modes (mesons). 

6.2. Plasma-Ball Production in Hadron-Hadron Collisions 

Consider the collision of two stable glueballs (with masses of order Kgap) at center of 
mass energies large compared to N^Agap. In the large A string theory dual these glueballs 
map to light modes such as gravitons (or other light flelds). Giddings has emphasized 
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that, at large A, black hole production saturates the inclusive cross section for graviton 
collisions at high enough energies (up to numbers of order unity); furthermore, this cross 
section also saturates the Froissart unitarity bound. We pause to review this argument 
(closely related to an argument given by Heisenberg in 1952 using pions instead of 

gravitons). An upper bound for the inclusive scattering cross section is easily estimated 
by assuming that gravitational forces dominate at high enough energies. Then, the force 
between two colliding particles separated by impact parameter b may be estimated to be 
of order '^'^f g-^gapb^ where we have used the fact that gravity is massive in the relevant 
backgrounds. Consequently, incident particles simply sail past each other when b is larger 
than a number of order \ii{E) / Agap (up to corrections that are subleading at large energies), 
and the inclusive cross section is bounded from above by 

a ~ ln\E)/Al,^. (6.1) 

However, one may independently use shock wave metrics and singularity theorems to 
demonstrate that the cross section for black hole formation is of the order of a in ( |6.1j ) 

m- 

As we have explained above, localized black holes map to localized lumps of gluon 
plasma in the dual field theory. It follows as a prediction of the AdS / CFT correspondence 
that, at large and large A, two glueballs shot at each other with an impact parameter 
smaller than (ln{E)/Agap) will coalesce into a lump of gluon plasma with a probability 
close to one when E ^ N'^Agap. This lump of plasma, produced at a rate that saturates 
the Froissart bound, will typically settle into a very long lived plasma-ball, which proceeds 
to slowly hadronize over a time scale of order A^ , in the manner described in the previous 
subsection. Note that even though the black hole forms with a size of order ln(i?) in the 
spatial directions of the field theory, it quickly expands to a size of order E^^^ (when the 
dual field theory has p spatial dimensions) for which it can be meta-stable. 

The reader may find the picture sketched in the last paragraph clashing with her 
QCD-trained intuition, which might lead her to expect the fast partonic constituents of 
the glueball to either pass right through each other or to undergo a small number of hard 
collisions rather than smoothly coalescing into a plasma-ball. Of course the utility of 
partonic ideas is questionable at large A (where partons interact strongly at all energies). 
Nevertheless, to the extent that this notion is valid, rapid gluon radiation ensures that the 
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parton distribution function is peaked at small values of x [T^\ so that the average parton 
energy is of order 

Eparton ~ E{—^^)^ (6.2) 

where E is the center of mass energy of the collision, ensuring that Eparton ^ ^gap at 
large A, so that glueballs simply do not contain fast partons. Thus, at large and large A, 
glueball-glueball collisions are conceptually similar to heavy ion collisions, with the large 
center of mass energy shared between a large number of constituents. 

At small A, as mentioned above, most of the energy of the glueballs is carried by a 
small number of partons, each of which carries a significant fraction (an energy of order 
E/ ln(i?)) of the energy of the glueball (at high energies the partons tend to have smaller 
values of x due to asymptotic freedom, but this is a logarithmic effect that does not affect 
our arguments). These very energetic partons interact weakly; their interactions will not 
form a plasma. As a consequence we expect a crossover in the dominant behavior of high 
energy scattering at A of order one. HI 

It would be interesting (and may be possible) to quantitatively verify the qualitative 
picture sketched in this subsection. 



7. Black Hole Physics from Gluon Plasmas 

As we have argued above, in confining gauge theories at large A, black holes in the IR 
of the background ( |1.1| ) are dual to plasma-balls. As discussed below, the decay process 
of these plasma-balls goes through ten dimensional Schwarzschild black holes. To our 
knowledge this is the first proposal for the dual of classically stable small black holes 
(black holes whose size is small compared to the curvature scale of the space that hosts 
them). As a consequence, we are able to use the setup of this paper to inquire how various 
mysterious phenomena involving black holes (for instance, the complete absorption of any 

Note that even at small A the glueball still contains a large number of small x partons whose 
strong interactions might create a lump of gluon plasma. However, this lump will not carry a 
finite fraction of the center of mass collision energy. One might think that even if only a small 
plasma-ball is formed, the fast partons would be bound to it by a string (since they carry a color 
charge) so they would eventually be pulled back into it, given that strings cannot break in the 
large N limit. However, since in order to form the plasma we need an energy at least of order A'"'^, 
the fast partons will also carry an energy of order A^^, allowing the relevant strings to snap even 
in the large N limit (since their length would be of order N"^). 
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incident object) manifest themselves from the dual viewpoint. The dual picture also allows 
us, in principle, to investigate how a' and gs corrections modify the classical properties of 
black holes. In this section we will present a very preliminary discussion of these extremely 
interesting questions. 

1.1. The Fate of Small Schwarzschild Black Holes 

As the pancake-shaped black holes (at large A) lose energy into radiation they shrink 
in size until their size approaches the curvature length scale of the space ( |1 . 1|) that hosts 



them; they then localize on the internal manifold (via a Gregory-Laflamme [|T7| , p!8[] type 
transition)ii, and then shrink further until their size is much smaller than the length 
scale of the background curvature. At this point the black holes resemble ten dimensional 
Schwarzschild black holesil. Thus, the subsequent evolution of these black holes is identical 
to that of black holes in flat ten dimensional space. As these black holes continue to lose 
energy by Hawking radiation, they eventually shrink to the string scale; the qualitative 
nature of their subsequent evolution is of great interest (see appendix D for a discussion). 

Of course the slow decay of warped black holes, described in the paragraph above, 
has a dual description as the decay of plasma-balls (by hadronization) in gauge theories at 
large A. The properties of small (string sized) black holes map to the properties of small 
plasma-balls, a reformulation that may be useful. 

7.2. Information Conservation in Hawking Radiation 

It has often been pointed out that the AdS/CFT correspondence ensures that black 
hole evaporation is a unitary process. Our identiflcation of speciflc gauge theory con- 
figurations (occurring in the decay process of plasma-balls) with 'small' ten dimensional 



Note that such a transition is not expected to occur at small A, and this may cause some 
differences between the details of the decay process at small A and at large A. Note also that this 
localization transition may not be smooth, so part of the energy of the black hole may be lost in 
the transition, but we expect the end-point to be a localized black hole carrying a finite fraction 
of the initial energy. See [24| for an initial study of the localization transition on AdSs x . 
Note that there is still some question whether the Gregory-Laflamme instability leads, in finite 



time, to an end state with localized horizon [25|. See |26] and references therein for a study of the 
Gregory-Lafiamme transition from a dual gauge viewpoint, in a similar context. 

Note that this is true for all known theories with large A, despite the very different UV 
physics of these theories. It would be interesting to understand this universality better from the 
field theory point of view. 
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Schwarzschild black holes makes this argument more specific. The production of a plasma 
ball in a glueball-glueball collision, and its subsequent decay via hadronization is clearly a 
unitary process; the end point of this process is only approximately thermal. 

Of course a full resolution of the information paradox requires the identification of 
the flaw in Hawking's argument that predicts a breakdown of unitarity. While the dual 
description of Hawking radiation is manifestly unitary, it has not yet proved possible to 
formulate Hawking's argument in gauge theory language in order to identify its flaw. This 
is an important problem that deserves attention. 

1.3. How Black is a Black Hole ? 

It is a striking feature of classical black holes in general relativity (the feature respon- 
sible for their name) that a particle squarely incident on the black hole is always absorbed, 
no matter how large its energy. 

In this paper we have identifled a black hole with a plasma-ball, a localized lump of 
gluon plasma at temperature T^. It might, at flrst, seem that a particle incident on the 
plasma-ball at an energy E much larger than would sail right through it (undergoing a 
series of grazing collisions that hardly affect it), in blatant contradiction with the expec- 
tations of black hole physics. We have already explained in subsection 6.2 that this naive 
expectation is incorrect. The only objects available to be hurled at the plasma-ball are 
glueballs. At large A glueballs may be thought of as consisting of a very large number of low 
energy partons; the average parton energy is Eparton in (|6.2|) . Note that Epart ^ ^gap even 
in the limit E ^ oo. As a consequence, glueballs may always be absorbed by plasma-balls, 
no matter how high the incidence energy. 

As we have explained in subsection 6.2, however, glueballs at small A are composed 
dominantly of a small number of high energy partons. We expect that these partons should 
simply blast through the plasma ball at high enough E (note that E may have to be of 
order in order to permit these partons to snap conflning strings that may bind them to 
the plasma-ball and to escape all the way to inflnity, since the snapping probability per 
unit time and unit length of the string is of order 1/N^). Thus, plasma-balls at small A 
(where string corrections are very important) do not share all the properties of their large 
A cousins (in particular they no longer efficiently absorb glueballs). 

The results of this subsection were obtained in collaboration with Nima Arkani-Hamed. 
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It is tempting to translate the last paragraph into a statement about a reduction in 
the absorptive properties of black holes in warped backgrounds whose curvatures are of the 
order of the string scale. Of course, this translation has an obvious pitfall; a glueball blast- 
ing through a plasma-ball may have a more mundane dual description, as a graviton that 
evades the corresponding black hole by going around it in the radial directionll (needless 
to say, all these notions are a little fuzzy when curvatures are stringy). Nonetheless, this 
suggestion is intriguing and deserves further investigation. 



8. Numerical Solutions for Domain Walls in some Specific Backgrounds 

8.1. Scherk-Schwarz Compactifications of Anti-De Sitter Space 

In this section we turn to a detailed study of the specific warped background 

ds^ = L'^a' ( e^" (-dt'^ + T2^(u)de^ + dwf) + —^—-duA , (8.1) 
V T2^{u) J 

where i = 1, ■ ■ ■ ,d — 1, 6 = 6 + 2ti and 

TAu) = l-[-{d + l)e-) . (8.2) 



This metric, known as the AdS soliton , is a solution to the d -f 2-dimensional Einstein 
equations with a cosmological constant 

and has a simple physical interpretation [§. It may roughly be thought of as a Scherk- 
Schwarz compactification of AdSd+2 on a circle. Indeed, at large w, Tr(w) ~ 1 and ( |8?1] ) 
reduces to AdSd+2 in Poincare-patch coordinates, with u as the radial scale coordinate, and 
with one of the spatial boundary coordinates, 9, compactified on a circle (the remaining 
boundary coordinates, Wi and t, remain non-compact). At smaller values of u, ( |8.1J ) 
deviates from the metric of periodically identified AdSd+2- In particular the 9 circle shrinks 
to zero at a finite value of u, smoothly cutting off the IR region of AdSd+2- 

In order to describe thermal physics about ( |8.1|) it will be convenient to switch to 



Euclidean space. Compactifying time r = r + /3 on the Euclidean continuation of (|8.1 



ds^ = L^a' f e^^ (dr^ + T2Md9^ + dwf) + —^duA , (8.4) 
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This was discussed in a somewhat different context in |27]. 
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we obtain the Euclidean configuration for a tliermal gas of gravitons (at temperature 



T = about dO]). The background 



ds^ = L^a' ( e'" (Tg{u)dT^ + dO^ + dwf) + -^dvA (8.5) 

yields a second smooth Euclidean manifold with the same asymptotics; upon continuing 
to Lorentzian space this solution has a horizon, so we identify it with the black brane at 
temperature T = 1/(3 about (^]l|) . 

It is immediately evident that for P = 2n the solutions (|8.4|) and (|8.5|) (the thermal gas 
and the black brane) are identical (they differ only by a labeling of circles) and so have the 
same free energy. In appendix E we compute the boundary stress tensor and free energy for 
these two solutions at every temperature, and demonstrate that the thermal gas has lower 
free energy for P > 2n, while the black brane solution dominates the thermodynamics 
at all higher temperatures. We also explicitly verify that the pressure of the black brane 
vanishes at the phase transition temperature Td = l/27r, as anticipated by the general 
arguments of section 4. 



Fluctuations around (|8.1| ) have a mass gap, and the phase transition temperature is of 
the same order as the mass gap (both are of order one in our conventions). It is also easy 
to check that the thickness of the brane at the phase transition temperature (the minimal 
value of u in ( ^.5| )) and the confinement scale of massive particles (the exponential decay 
in the radial direction in solutions of = in the background ( |8.1| )) are both of order 
one. 

Therefore, the background ( |8.1|) is a particularly simple warped background of the 
form described in subsection 1.2; its simplicity stems partly from the fact that it is a solu- 
tion to Einstein's equations with a cosmological constant but no additional matter fields. 



Moreover, the background (|8.1|) also appears as a component of string compactifications 
that are dual to interesting field theories, as we now review. 

Some backgrounds of the form ( |8.1| ) xM are solutions of string theory (with appro- 
priate fluxes) that are holographically dual to d + 1 dimensional conformal field theories 
with one spatial dimension compactified on a circle with anti-periodic (Scherk-Schwarz) 
boundary conditions for the fermions. When d = 3 the field theory dual to the simplest 
such compact ificat ion (that with M = S^; different choices of M lead to different field 
theories) is AT = 4 supersymmetric Yang Mills theory compactified on a Scherk-Schwarz 
spatial circle 0. At low energies (and for small A) this theory reduces to pure Yang Mills 
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theory in 2 + 1 dimensions, a confining gauge theory, making contact with the gauge theory 
discussions in this paper. When d = 4 the field theories dual to such backgrounds are five 
dimensional superconformal theories compactified on a circle with anti-periodic boundary 
conditions for the fermionic fields. Such theories are effectively 'confining' in four dimen- 
sions, and in some cases the corresponding five dimensional conformal field theories may 
be viewed as the strong coupling limit of a five dimensional gauge theory [P^ . 

For all these reasons, in the next subsection we will turn to a detailed study of the 
background ( |8.1| ). In particular we obtain a numerical construction of the domain wall 
that interpolates between ( pl^ ) and ( |8.5|) at /3 = 27r. 

8.2. The Domain Wall: Asymptotic Behaviors 

As discussed in §4, the black hole solutions we are interested in should interpolate 
between the (Lorentzian versions of the) confining vacuum ( |8.4| ) on their outside and the 
black brane solution ( p.5|) at the deconfinement temperature on the inside. In the limit 
that these black hole solutions are very large, the solution in the neighborhood of their 
edge should reduce to a domain wall (in one of the field theory directions) interpolating 
between these two solutions. 

Since the metrics must depend explicitly on two coordinates (the bulk radial coordi- 
nate and the field theory direction normal to the wall), it is unclear whether any analytic 
solutions may be found. The elegant Weyl metric, whilst extended to higher dimension, 
has not been extended to include a cosmological constant [p0|,|31|j3^ . Hence we will use nu- 



merical methods to construct the metric, following those employed previously in ||33| , p4| , p5 
(see also [^). Note that whilst it seems possible that the domain wall system we solve 
here may one day find an analytic solution, similar equations with more bulk matter, such 
as that found in general confining dual backgrounds , would likely not be solvable, 

whereas the numerical methods we set-up here should easily extend to these cases. Fur- 
thermore, the finite size black hole problem is likely to remain analytically intractable.0 
We focus our calculation on the domain wall solution rather than on finite size solutions 



There has been interesting analytic progress in lower dimensional (d = 2) cases with just a 
negative cosmological constant and the bulk ending at a thin brane |]37|j3^] . This progress follows 



from the AdS C-metric [39|, which unfortunately is not known in more than 4 dimensions. We 
also emphasize that everything we say applies only to warped backgrounds with a mass gap (and 
with no four dimensional massless gravity). In particular the analysis of this paper may not apply 
to 'cut off' AdS geometries with a massless 'radion' field. 
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simply because the equations and solution are most elegant in this case. Moving to finite 
size is in principle a straight forward extension of the methods we use, which we leave for 
future work.i^ 

In this subsection we give an overview of our numerical construction of this domain 
wall. The formal construction, detailing the equations to be solved, the boundary condi- 
tions and the data, is presented in appendix F. Numerical details are given in appendix 
G. We are mostly interested in the cases of (i = 3, 4 but we present our formulae for ar- 
bitrary dimension. We will consider the problem in the Euclidean setting, although we 
stress that the equations, and their boundary conditions, are totally independent of the 
spacetime signature. The static Lorentzian solution is obtained by trivial continuation of 
the Euclidean one. 

We will search for a solution to ( ^.3|) that preserves rotational and translational sym- 
metry in {d — 2) of the {d — 1) spatial Wi directions, which we call (a = 1, ■ ■ ■ , ((i — 2)), 
and translational symmetry in the 9 and r directions. The background is also required to 
respect refiection symmetry in r and 9. It follows that it is locally possible to choose the 
metric for our background to take the form 

ds^ = A^dr"^ + B^d9^ + e^^drl + e^^ {dx^ + dy"^) , (8.6) 

where the metric functions depend on the coordinates x, y which represent some combina- 
tions of the field theory direction on which the domain wall is localized and of the radial 
coordinate in the bulk space (we choose specific combinations such that the metric takes 
the form (P^)). We take the coordinate periodicity of the r, 9 circles to be /3, (3 respectively, 
although we note that in our asymptotic AdS space-time (or, from the perspective of the 
{d + 1) -dimensional field theory, due to conformal invariance in the UV) the actual value 
of /9 or /3 is irrelevant. It is only the ratio of the proper size of the Euclidean time circle 
in the UV to the proper size of the space circle in the UV that is physical. As discussed 
above, these are equal exactly at the deconfinement temperature and we only expect to 
find a Euclidean (and hence static Lorentzian) solution in this case. 

Note that the ansatz (|8.6[) leaves a residual diffeomorphism freedom of performing 
conformal coordinate transformations of the x,y plane which preserve the form of ( |8.6|) . 

This will introduce new terms in the Einstein equations as planar symmetry in the field 
theory directions is replaced by spherical symmetry, and similarly the boundary conditions will 
then include the origin of these spherical-polar coordinates. However the same methods will apply, 
and in particular the metric will still depend non-trivially on only two coordinates. 
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Figure 4: Domain of problem in x,y plane. 

Remembering that (non-singular) two dimensional conformal transformations can map 
any region into any other region (preserving angles at corners) , we will use this freedom to 
position the boundaries of the domain of our coordinates x and y as in figure 4. We will 
take the IR boundaries to be at x = and y = 0, and the UV boundary will be defined by 
the vanishing of a function /(x, y). For large x we take the function / to vanish at y = c, 
for some constant c. For large y, similarly we take / to vanish at x = c. Furthermore we 
will choose / so that the function is invariant under the symmetry /(x,y) = /(y,x) 
(which will yield a symmetry of our solution, to be discussed below). We regard the 
locus of the zero of /, and thus the particular value chosen for c, to be fixed throughout 
the following discussion. Note that the particular value of c chosen is irrelevant, and may 
be changed simply by performing a global scaling {x,y) (Ax, Ay), which rescales the 
UV boundary position and the metric functions A, B, e'~^ together with /3, (3 but does not 
physically change the solution. 
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The metric may always be written in the form ( p.6| ) locaUy. We assume that there 
exists a solution for which the metric functions A,B,C,D in ( |8.6D are finite everywhere 
in the coordinate domain of the problem, except at the conformal boundary in the UV 
where / has a first order zero and A, S, e*-", ~ 1//. It is certainly not obvious that such 
a solution exists, but we will show that it does by constructing it numerically. Note that 
once we have made the choice of coordinate boundaries above, we have completely fixed the 
residual conformal coordinate freedom (up to pathological transformations which diverge 
at large x,y and thus are projected out under our finite metric component assumption). 

Let us discuss the boundary conditions in more detail. The bulk space-time has a 
conformal boundary at f{x,y) = 0, and the region near / = looks like anti-de Sitter 
space near its conformal boundary. Thus the conformal boundary metric is (conformally) 
flat. We require the bulk to close smoothly in the IR, with either the r or 9 circles 
shrinking to zero size. As indicated in figure 4 we choose the time (r) circle to shrink for 
X > 0,y = and the space (9) circle for y > 0, x = 0. At the origin x = y = 0, both circles 
smoothly vanish, as discussed in more detail below. At large x and at large y the geometry 
should tend to one of the homogeneous solutions (|8.4|) and ( |8.5|) with the appropriate circle 
shrinking, as detailed in the previous subsection. 



The residual conformal transformations that preserve the form of ( p.6| ) and maintain 
the finiteness of the metric functions must preserve angles at corners of boundaries of the 
domain. Hence we have made a crucial choice above, namely that the angle between the 
X and y axes is a right angle. Let us now justify this choice. If we focus on a small enough 
region near the origin in x and y we may ignore the bulk cosmological constant. The local 
geometry here must be flat space, since we require the Euclidean geometry to be smooth. 
Hence, the metric must take the form 



A^dr^ + B^d9^ + {dA^ + dB^) + [e^^dr^ + 0{A^, B^) (8.7) 



to leading order in A, S around A = S = 0, where to leading order C is a finite constant, 
and r, 9 are rescalings of r, 9 with coordinate period 27r, such that this Euclidean manifold 
is smooth. This is simply [d + 2) -dimensional fiat space, written as a product of four 
dimensional fiat space, in double polar coordinates, and the fiat {d — 2) -dimensional space 
parameterized by the directions. In ( p.7| ) the lines defined by the shrinking circles, A = 
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and B = 0, are indeed at right angles, and the metric components are finite, justifying our 
choice for the IR boundary.!! 

The analytic details of the construction are found in appendix F. There we find that 
the boundary conditions have no remaining data once we have specified our IR and UV 
boundary locations (and fixed the trivial rescalings of the r, 9, coordinates). The equa- 
tions and boundary conditions in the IR and UV (having chosen appropriately the UV 
boundary location as described above) are invariant under the symmetry (x, y) ^ [y, x) 
providing we interchange A and S, and the asymptotic solutions ( |8.4|) and ( ^?5|) we expect 
to find at large x and at large y are also exchanged by the symmetry. Hence we expect 
that the solution for the domain wall will possess this 'Z2 symmetry so that 



For a smooth Euclidean geometry this implies /3 = /3, so the time and spatial circles have 
precisely the same proper size in the UV, and so the Euclidean temperature is indeed 
exactly at the deconfinement temperature, as discussed in the previous subsection. Al- 
ternatively, in the Lorentzian setting, the surface gravity of the horizon, as measured by 
an IR observer far from the domain wall in the confining region, is directly related to the 
circle size in that region. 

Having fixed the coordinate boundary locations (and fixed the trivial rescalings of the 
T,9,ra coordinates), we obtain a unique solution for the metric functions A,B,C,D in 
both (i = 3, 4 with the above symmetry. We may then determine the periodicity /3 of the 
two circles from the proper gradient of the vanishing of the dr"^ metric component in the 
asymptotic black brane region (or equivalently, the vanishing of the d9^ component in the 
confining region) . We see in appendix F that the Einstein equations indeed guarantee that 

Define z = A -\- iB, and consider a singular conformal transformation that takes the form 
w = near z = 0, with w = A + iB. Such a transformation with p 7^ 1 will change the angle 
between the two axes. The new metric will contain the terms \wY'^'^{dA^ + dB'^) from the Jacobian 
factor of the transformation, so e^^ = For p 7^ 1 this means \D\ is infinite at z = w = Q. 

Hence if we require a metric of the form ( |8.6| ) with all the metric functions being finite, we must 
have a right-angle between the axes defined by the shrinking of the Euclidean time circle, and 
that of the space circle. 



A{x,y) 
C{x,y) 
D{x,y) 



B{y,x), 
C{y,x), 
D{y,x). 



(8 



8) 
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Figure 5: Surface plot of / in the x,y plane (in fact we plot (— /) for 
visual convenience) for d = A. The zero of / gives the position of the UV 
boundary. / is not defined over our entire coordinate domain, but rather in 
a finite region near the UV boundary. 

the circle coordinate periods required to smoothly close the appropriate shrinking circle in 
the homogeneous regions indeed smoothly close the circles everywhere in the solution, in 
close analogy with the zeroth law of black hole mechanics. 

In figures 5 and 6 we present the form of the metric functions we compute numerically 
for the case d = 4. For d = 3 their qualitative form is very similar. We work in units 
where L^a' = 1. In figure 5 we show the function /, whose zero gives the position of the 
UV boundary. Note that this function is not defined over the entire domain, as it is only 
required in a neighborhood of the UV boundary. The details concerning the functional 
form of / are discussed at length in appendix F. Then, in figure 6 we show surface plots of 
the metric functions A, S, e*^, as functions of x, y. We see that these behave smoothly 
everywhere, except near the UV boundary where they diverge as ~ 1/f (where / has a 
first order zero). A goes to zero at y = 0, and hence the time circle shrinks there. By 
symmetry, B goes to zero at a; = where the space circle shrinks. We clearly see that at 
large x (and by symmetry large y) the solution does indeed quickly become homogeneous 
as expected. 

We compute the boundary stress tensor for both the d = 3 and d = 4 solutions in figure 
7. The horizontal axis is the proper distance along the boundary (up to a global scaling 
explained below), which we identify with a coordinate g in the field theory perpendicular 
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Figure 6: Surface plots of A, B, e , e in the x,y plane for d = 4. The 
corresponding plots for d = 3 appear very similar. 

to the domain wall. The vertical axis gives the various independent components of the 
traceless stress tensor. A key point is that Tgg is constant (in fact it is zero), since the 
stress tensor conservation yields dgTgg = (we discuss this in detail in appendix F). 
As discussed in appendix B.l, the surface tension of the domain wall is given by 

/c» 
dgTr^r^ (8.9) 
-oo 

(no sum over a). From our numerical results depicted in figure 7 we find that the domain 
wall tension E, measured in units (as in the figure) of the energy density pc of the deconfined 
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Figure 7: The components of the boundary stress tensor (in units where 
T-rr — > 1 in the deconfined phase) as a function of a coordinate g perpendicular 
to the domain waU (negative is confined region, positive is deconfined region), 
for (i = 3 (left) and d = 4 (right): red (top) = T^t-, orange (bottom) = Tgg, 
green (top middle) — T^^ra^ blue (bottom middle) = Tgg. 

phase at the deconfinement temperature divided by the deconfinement temperature 
(which was used to set the scale of the g coordinate in figure 7), is given by E = 2.0 x Pc/T^ 
for d = 3 and by E = 1.7 x pc/Td for d = 4. We estimate the systematic error in these 
numbers to be no more than ten percent. 



K 




1 4 



Figure 8: Surface plot of the scalar curvature invariant K = {Rfj,i,jpR^'^'^P)^/^ 
X, y plane for d = 4. 
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Finally, in figure 8 we plot the scalar curvature invariant (i^^j^^pi?'*^''''')^/^ for the d = A 
solution. We see that near the UV boundary it tends to a constant, namely that for AdS. 
In the homogeneous regions it is a little larger in the IR, reflecting the fact the solution 
is really Euclidean AdS-Schwarzschild, with greater curvature localized where the circle 
shrinks. In the domain wall region, we see the curvature is increased, but it is everywhere 
smooth as we expect. 
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Appendix A. Thermodynamics of Large N Gauge Theories 

In this appendix we review the thermodynamic properties of large N SU (N) gauge 
theories that undergo first or second order deconfinement phase transitions. We discuss 
exphcitly theories having only fields in the adjoint representation (such as pure Yang-Mills 
theory and its supersymmetric generalizations), but the discussion does not change if we 
add a finite number of fiavors. 

A.l. Local Stability of Homogeneous Configurations 

We begin by reviewing the thermodynamic criteria for the stability of a homogeneous 
system. 

Consider any system (like a gauge theory) described by a local field theory. In the 
micro-canonical ensemble the system is characterized by the energy density p = E/V, and 
the thermodynamical behavior is determined by the entropy density f{p),S = Vf{p). The 
effective temperature is given by T = l/{dS/dE) — l/f'{p). Here we have assumed that 
the system is homogeneous, but of course the homogeneous configuration does not always 
maximize the entropy. To see this, divide the system into two pieces, of volume aV and 
(1 — a)V, respectively. Let the energy in the two parts be {a + 5a)E and {1 — a — 5a)E, 
respectively. Sa = corresponds to the homogeneous system; at small non-zero Sa the 
entropy of this system, above that of the homogeneous configuration, is given by 

5S _ 5a'p' d'f 

V 2a(l - a) dp^ ■ ^ ^ 

Note that 5S has the same sign as d'^f/dp'^. It follows that the homogeneous phase is 
locally unstable when d'^f/dp'^ is positive. Since df /dp is the inverse temperature of the 
system, d'^f/dp'^ has an opposite sign from the specific heat of the system. Consequently 
we have merely rederived the well-known fact that homogeneous systems with negative 
specific heat are unstable. 

Note also that the pressure of a system is given hj p — —dE/dV\fixed s- Using 
dS = dVfip) + f'{p){dE - pdV), we find that 

p{p) = -^^jr^ = -Ffree{p), (A.2) 

where Ffree is the free energy density. The velocity of sound squared is given by 

2 ^dp^ zlll (A 3) 

'^sound (/')^ " V / 

It follows that small pressure waves in this homogeneous system are tachyonic (so that the 
system is dynamically unstable to small perturbations) if and only if /" > 0, reiterating 
the conclusion of this subsection. 
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A. 2. Density of States in Confining Gauge Theories 

We now specialize to the study of a confining large Yang-Mills theory that undergoes 
a first order deconfining transition, and review the behavior of the function /(p) in this 
class of theories (which may include the models of [l8|JT^,p|,[T0(| at various values of A) . Two 
interesting mass scales in such theories are Th ~ Ms, the Hagedorn temperature of the 
low energy confining large N gauge theory (related to the confining string tension M^), 
and the mass gap Agap. At large A, Th / Agap is proportional to a positive power of A, while 
for small A this ratio is a number of order one. 

We first study /(p) at energy densities which are 0{N^) (rather than 0{N'^)). When 
p <^ Mg the system is well-approximated by a gas of weakly interacting glueballs, a system 
that has positive specific heatH. When the energy density reaches the string scale, p ^ T^, 
the system is well-described by a gas of weakly interacting strings at approximately the 
Hagedorn temperature Th, f{p) oc p/Tn (this is a consequence of the exponential or 
Hagedorn growth in the high-energy density of states e^/^"). Note that cP f /dp^ vanishes, 
indicating marginal stability, when p ^ Tfj 0. 

Next, we turn to energy densities of order A^^. We first consider the very high energy 
limit. At very high energy densities the system is governed by its high-energy 'fixed 
point'. This fixed point could, for instance, be a free theory with a number of particle 
species of order A^^, or a conformal field theory in four dimensions. In either of these 
cases f{p) oc N^^'^p^/'^ for p ^ N'^A^^^ Ell. In this regime the effective temperature grows 
indefinitely with increasing density and the system has positive specific heat. 

We now turn to intermediate energy densities of order N'^. When the deconfinement 
transition is of first order, it must occur at a temperature T^ lower than Th (otherwise 
the deconfinement transition would be a second order transition at Th] clearly it does not 
make sense for the confining phase to exist at temperatures larger than Th)- Thus, there 



In theories with a single mass scale, Ms is of order Agap and this energy range is pretty 
boring. However, at large A, Ms is parametrically separated from Agap and this energy range is 
more interesting - typically the density of low mass glueball states grows as a power of the mass 
(due to Kaluza-Klein modes in extra dimensions), leading to f(p) oc with 3/4 < a < 1. 

More precisely, the behavior in this regime depends on the sign of the corrections to the 
Hagedorn behavior, f(p) ~ p/Th +Qln(p), but in any case the second derivative is very small for 
large p. 

Some of the large A theories that we are interested in actually involve more exotic UV 'fixed 
points', which have a different form of /(p) but the same qualitative behavior. 
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must exist an intermediate regime between the Hagedorn regime and the high-energy- 
density regime, in which the effective temperature decreases to a value below T^. Such 
a phase obviously has negative specific heat and is unstable S As p increases in this 
phase, we expect that d? f / dp^ decreases monotonically, passing through zero at some 
p = po of order A^^A^^^, and eventually joining with the negative values of d'^ f /dp^ in the 
high-energy-density phase. 

-Log(S'(E)) 
= Log(T) 



Log(T_d) 
Log(T_0) 



Log(EA') Log((E/V LO) Log((E/V)_l ) 

Figure 9: The (logarithm of) the effective temperature T = l/{dS/dE) as 
a function of the (logarithm of) the energy density in large A^ confining gauge 
theories which undergo a first order deconfinement transition. 




The behavior described above is depicted in figure 9, which plots the logarithm of the 
effective temperature as a function of the logarithm of the energyH. Stable homogeneous 
phases correspond to lines with a positive slope in this figure. 

The generic behavior of systems exhibiting a second order phase transition is quite 
similar, except that the unstable intermediate phase does not exist. In such systems the 
Hagedorn phase joins smoothly with the high-energy-density phase and there is a second 
order phase transition precisely at the Hagedorn temperature 



"^^ In many gravitational duals at large A, this phase is dominated by Schwarzschild-like black 
holes in ten dimensions, for which /(p) oc N^'^^'' p^^'' . 
"^^ A similar behavior appears also in the gravitational dual of (p + l)-dimensional conformal 



field theories compactified on [iO,|]; see [41| for a recent discussion. 
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A. 3. Phase Separation in Confining Large N gauge theories 

Consider the behavior of a large gauge theory which has a first order deconfinement 
transition, at an energy density p which is of order A^^ but smaUer than pq. As discussed 
above, the homogeneous phase is unstable in this range of densities, so we expect the 
system to decay into separate regions in two different phases, one with a density smaller 
than p and the other with a density larger than p. We expect the two phases to be the 
stable phases in the discussion above, which are the low-energy-density confined phase, at 
p < Tfj, and the high-energy-density deconfined phase, at p > N'^Ag^p. In the first phase 
the energy density, entropy density and free energy density are all 0{1). We will find it 
convenient to normalize the extensive quantities (such as E and 5") by dividing them by a 
factor of A^^ as well as by the volume. All normalized extensive quantities vanish in the 
first phase. 

The second phase in this non-homogeneous mixture must have an energy density 
larger than po in order to be stable. We now determine the properties of this phase. Let 
a fraction a of the net volume of our system be occupied by this phase; since the energy 
density of the confined phase is negligible, this means that the energy density in the second 
phase should be p/a. Now, a may be determined by maximizing the entropy of the system, 
which is done by by maximizing af{p/a). This is the same as maximizing S{p')/p', where 
p' is the energy density of the second phase. The energy density p' = p/a of the second 
phase is determined by this maximization, which gives the equation 

f'ip') = (A.4) 
P 

It is not difficult to see that equation (|A.4| ) has a single solution pi, and that pi > po. 
Equation ( |A.4| ) may be physically interpreted in several ways. First, recall that the 
free energy density of our system is given by Ffree{p) = {E — TS)/V = p — {f{p)/f'{p))- 



The argument goes as follows. Recall that, for p < po (but still of order A'^'^), /' is an 
increasing function. It follows that, for p < po, f'{p) is more than its average value {f{p)/p) 
between and p. Consequently, (|A.4| ) has no solutions for p < po. Now, for p in the high-energy- 



density phase, the right-hand side of (A.4) is larger than its left-hand side. Consequently, ( |A.4D 
must have at least one solution for p > po- Let the smallest such solution be pi. Using the fact 
that /' decreases for p > pi, it is easy to argue that no p > pi can obey ( [A.4 ). It follows that pi 
is the unique solution to ( [A. 41) . Note that we assume here that there is only one point where the 
second derivative of / vanishes (as in the discussion of the previous subsection); this is true in all 
known cases, though more complicated possibilities cannot be ruled out. 
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Comparing with ( |A.4|) , we conclude that our system has positive free energy density for 
p < pi but negative free energy density for p > pi; the free energy density vanishes at 
p = pi. As the (normahzed) free energy density also vanishes in the first phase, the free 
energy densities of the two phases are equal when the deconfined phase has p = pi, so we 
conclude that pi is the energy density of the gluon plasma just above the deconfinement 
temperature. Second, recalling that the pressure is simply equal to minus the free energy 
per unit volume ( |A.2|) , it follows that p{p) is negative for p < pi but is positive when 
p > Pi, and the pressure vanishes precisely aX p = pi. 

The bubbles of gluon plasma at density pi which appear in this discussion resemble 
the localized bubbles discussed in section 2; the main difference between them is that the 
bubbles of section 2 had the vacuum outside them rather than the finite-density confined 
phase, but this difference is negligible in the large limit. 



Appendix B. Properties of Plasma Dynamics 

B.l. The Surface Tension of the Plasma-Ball 

The force balance equation (|2.1| ) may be derived from stress energy conservation 
applied to the plasma fluid in polar coordinates. In p spatial dimensions the equation 
V^T^r = 0, for static configurations, reduces to 

drTrr + \ [pTrr - 9^'%, ) = 0, (B.l) 

where i,j are summed over all spatial coordinates. When the plasma fluid is isotropic 
pTrr = g^'^Tij, and constant Trr solves (|B . 1|) . In the plasma-ball configuration, however, 
the stress tensor will not be isotropic near the boundary of the bubble. Integrating ( p . 1| ) 
we find 

dr 

Trr{^) - TrM = P{0) - P{oo) = - / — {pT^r ' Q^'T,^) , (B.2) 

Jo ^ 

where -P(O) and P(oo) are, respectively, the pressure in the interior and exterior of the 
bubble. In the plasma-ball configuration the integral on the right-hand side of ( [B.2D receives 
contributions only from the neighborhood of the boundary of the bubble. Comparing 
with (|2.1| ) we conclude that the domain wall tension (for bubbles that are much larger 
than the width of the domain wall, which is expected to be of order l/A^^p, so that r is 
approximately constant in the region that contributes to the integral in (|B.2|) ) is given by 

E = / dr [g'^T.j - pT^r) ■ (B.3) 

38 



In the limit of an infinite bubble whose domain wall is transverse to the x direction, stress 
energy conservation implies that T^x is a constant (which vanishes since it is zero in the 
confined vacuum), and the expression for the domain wall tension reduces to 

■1 /"OO POD 

S = — - / dxg'^T.j = / dxTyy, (B.4) 
P ^ J — 00 J — 00 

where y is any spatial dimension along the surface of the domain wall. 

In this paper we have assumed that the effective dynamical surface tension between 
the confined and deconfined phases is positive near the phase transition temperature. This 
is true in the specific example analyzed in section 8. However, it may turn out that 
the effective surface tension is negative in some large gauge theories - at least we are 
unaware of an argument that rules out this possibility. Large plasma-balls will not be 
stable in such theories - they will be unstable to fragmentation into small plasma-balls 
(whose size is presumably of order the inverse mass gap). However note that, at least in 
gauge theories at large A, small plasma-balls (with R <^ 1/ Agap) are always stable as they 
map to Schwarzschild black holes in the dual gravitational theory. 

B.2. Dynamical Evolution of Lumps of Plasma 

Consider a lump of gluon plasma with large volume VA^^p ^ 1, at an average energy 
density larger than the critical density. The initial pressure of such a bubble is positive, 
driving the bubble to expand. As the expansion proceeds, the pressure and temperature 
gradients inside the bubble tend to even out. In large systems that undergo first order 
transitions we have mentioned in section 2 that the plasma-ball is an attractive fixed point 
in the space of dissipative fiuid fiows, and so it will be the equilibrium configuration if we 
start in its domain of attraction. We believe that the only other stable configurations of 
plasma fiuid are disconnected plasma-balls moving away from each other. It follows that, 
within fluid dynamics, hot lumps of gluon plasma stabilize into a collection of plasma 
baUsJl 

We will argue below that we do not expect stable plasma-fluid configurations to exist 
in systems that undergo second order deconfinement phase transitions. We expect the 
bubbles of plasma in such systems to expand through the deconfinement transition and to 
rapidly hadronize. 

Note that, under some circumstances, it is rather natural for a single lump of gluon plasma 
to break up into many components. This could happen if, for example, during the process of 
evolution the local energy density of some part of the fluid were to become low enough such that 
the local speed of sound squared is negative (it should go below po, see appendix A), triggering 
an instability that results in phase separation. 
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B.3. Small Plasma-Balls 

As discussed in section 3, plasma-balls shrink as they slowly lose energy by radiating 
glueballs. According to ( |2.1| ) the plasma-ball temperature diverges as its radius goes to 
zero. However, (|2.1D does not apply when R < 1/Agap, since the surface tension term 



does not capture all surface terms for small bubbles (contributions to the energy involving 
the curvature of the surface and inhomogeneities of the plasma are as important or more 
important). A more detailed microscopic analysis is needed to determine the properties 
of small plasma-balls. The strongly coupled nature of the gluon plasma may make this 
analysis rather difficult; however it could be worth the effort, as small plasma-balls encode 
very interesting information about black hole physics, as we have argued in §6 and as 
discussed further in appendix D. 

B.4- Plasma Lumps in Theories with Second Order Deconfinement Transitions 

Stable plasma-balls do not exist in theories that undergo second order phase transi- 
tions, at least when E is positive. While plasma-ball-like configurations of finite radius 
appear to exist as stable solutions of fiuid mechanics (using the arguments of section 2), 
their surface temperature is higher than the Hagedorn temperature, so they cannot be in 
an approximate equilibrium in a large gauge theoryll. Note also that in this case asymp- 
totically large plasma-ball-like configurations are unstable even within fiuid mechanics, as 
the specific heat diverges (and so the speed of sound vanishes) at the phase transition point 
in any second order transition. 

B.5. Plasma-Balls in the Real World ? 

A large N version of QCD may be obtained by coupling SU{N) pure Yang-Mills theory 
to (say) Nf = 2 or Nf = 3 light quarks (fermions in the fundamental representation). As 
we have remarked above, there is good evidence that this theory undergoes a first order 
deconfinement transition and so possesses plasma-balls at large enough A^. If the phase 
transition stays of first order as we decrease A, it seems possible that these meta-stable 
configurations will continue to exist (albeit with decreasing finite lifetimes). 



Recall that a hot ball above the Hagedorn temperature immersed into a large gauge theory 
cools instantaneously by radiating glueballs at a rate which diverges in the large limit. 

"^^ Note, however, that in QCD-like theories with Nf of order A^, plasma-balls would not have 
long lifetimes even if the phase transition remains of first order, because the number of different 
meson species that they can emit scales as . See for a recent discussion of the deconfinement 
transition for large ~ Nf. 
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The expected phase diagram of this theory at = 3, i.e. of real world QCD, may 
be found, for instance, in figure 1 of |^^. The nature of the deconfinement 'transition' in 
this diagram depends crucially on the strength of the chemical potential that couples to 
baryon number. When the chemical potential is small enough, QCD undergoes a smooth 
crossover between confining and deconfined behavior; this crossover becomes a strong first 
order phase transition at large chemical potential. As a consequence, QCD could possess 
meta-stable plasma-ball-like configurations at large enough baryon density, if the phase 



boundary surface tension is positive (see |44,45| for closely related discussions). However, 
meta-stable plasma-ball-like configurations will not exist at low baryon density, for instance 
at the baryon densities attained in the central rapidity region of the fireball created at 
RHIC, see for reviews. 

The rapidity tails of the RHIC fireball, as well as suitable regions of the fireball (if 
one was produced) in the SPS experiment, or plasmas produced in similar experiments 
where baryons are heatedH may have chemical potentials large enough to be in the first 
order regime. It would be interesting to analyze the relevant data from these experiments 
for evidence of formation of meta-stable plasma-ball-like configurations, that then decay 
into hadrons at approximately the phase transition temperature (see [ ^|50| and references 
therein for discussions along these lines). 



Appendix C. Counting Powers of at Lowest Order in Perturbation Theory 

In this brief appendix we illustrate the arguments of section 3 by explicitly displaying 
and computing the N dependence of sample diagrams that contribute to the gluon mean 
free path and to the glueball production rate at lowest order in perturbation theory. 

C.l. k = : Mean Free Path 

We check the result of §3.2 at weak coupling by studying the leading gluon-gluon 
scattering graphs shown in figure 10. First consider the graph in figure 10(A) for scattering 
an gluon and an (m, n) gluon, and let m ^ j and n ^ i. This diagram is proportional 
to g'^N^ = X/N, leading to a collision cross section proportional to /N"^ and so a mean 
free collision time that is independent of (the mean free time is obtained by multiplying 
the collision cross-section by the density of potential collision targets, whose number scales 
like N"^). 
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See figure 1 of |4S] for tlie location of various experiments on the QCD phase diagram. 
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A 



B 



Figure 10: Tree diagrams for gluon-gluon scattering in double-line notation. 

In the special case that (m, n) = {j, k) for some /c, the tree-level diagram in figure 10(B) 
also contributes. The cross section from this diagram is oc \^ /N, where the additional 
factor of N results from the sum over the q color index in the final state. This graph too 
leads to a mean free collision time that is independent of N as, in this case, the density of 
potential collision targets scales like N. 

C.2. k = 1 : Rate of Glueball Production 



Figure 11: Tree diagrams for glueball creation by scattering gluons. 

Consider glueball production via the gluon-gluon collision depicted on the left of figure 
11. The vertex in this graph contributes Qym^ index loop contributes a factor of A^, 
and the overlap of the two-gluon state with the glueball wavefunction is of order 1/N 
(this is the 1/N in the normalization of the glueball creation operator), resulting in a 
collision cross section of order l/N"^. Multiplying this by the total number of gluons in the 
plasma-ball, we find a glueball production rate independent of N. 

The same Feynman diagram may be given a different interpretation upon choosing the 
time slicing on the right of figure 11. This diagram describes an energetic gluon escaping 
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out of the surface of the plasma, until the string that attaches it to the rest of the plasma 
snaps by creation of a gluon-anti-gluon pair, allowing the original energetic gluon and its 
anti-gluon partner to escape away as a glueball. The analysis of the previous paragraph 
carries through immediately for this process, leading, once again, to a glueball creation 
rate independent of A^. 



Appendix D. The Final Decay of Small Plasma-Balls 

After the black hole/plasma-ball at large A shrinks such that it resembles a ten di- 
mensional Schwarzschild black hole, it continues to shrink as it evaporates until its proper 
Hawking temperature is of the order of the Hagedorn temperature Th (so that the proper 
size of the black hole is of order the string scale). At this point the black hole is believed 
to resemble a highly excited string state ||5^,^ ^ and there are several possibilities for 
the subsequent evolution. The black hole could 

1. Develop a tachyon (due to stringy a' corrections) at a temperature smaller than Th 
as a result of which the black hole decays classically, dominantly into low-lying string 
states (glueballs). 

2. Reach the Hagedorn temperature when its mass is of order 1/ g1 ~ A^'^, at which point 
it explodes into a gas of dominantly highly excited string states (or transmutes into a 
single excited long string, depending on details, see for instance section 2 of [0)^. 

3. Continue to decay thermally (with temperatures less than Th) until its energy is of 
order one. 

Which of the options 1-3 above is actually attained is a question of obvious interest 
in the study of black hole physics. This question could, in principle, be answered by an 
investigation of classical string theory, but this has not yet proved possible. Our identifi- 
cation of black holes with plasma-balls allows us to map this question to the properties of 
small plasma-balls, a reformulation that could lead to progress. 



"^^ See 1 41] for an intriguing discussion of a possible manifestation of the black- hole excited 
string transition in TV = 4 Yang Mills on . 

Note that the excited string states are stable in the large N limit, though for finite A'' they 
would eventually decay into lighter glueballs. 
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Appendix E. The Boundary Stress Tensor for Compactifications of AdSd+2 

Asymptotically AdSd+2 backgrounds possess a conserved d + 1 dimensional stress 
tensor |[5^ , |5^ , |55| , ^|57| . The stress tensor associated with ( |8.1D is (see equations 4 and 25 
in l5l 



Q,e-(c'-i)i' 

T^u = ©At!^ - - {d- l)-f^^ = (S^r, + {d + l)df,eSu9) , (E.l) 

where the indices /i and u run over i,t,9, 7^,^ is the restriction of (|8.1|) to a constant u 
slice, = —h }- — '^-hi^ and a = 1) ^'^^^^ , The boundary field theory stress tensor 
associated with ( pT|) is obtained by multiplying ( [E.l| ) by e^*^"^)", so that 

ycon/^ned vacuum ^ | ^ ^ \)b ^eKo) ■ (E.2) 

The boundary stress tensor for the black brane solution ( |8.5| ) is obtained similarly, giving, 

-^black brane /'o„^^^'^+l ^ 

2 



ybZacfc brane ^ (27rT)'^+^ - (-5^. + (d + l)d^tS.t) • (E.3) 



The subtraction we used to define the stress tensor in ( [iii.lj ) is the standard subtraction 
for the l)-dimensional field theory, which gives a vanishing stress tensor in the vacuum 
of this theory From the point of view of the d-dimensional field theory which we obtain 
after compactifying on the circle, it is more natural to use a different subtraction in which 
the stress tensor in the vacuum (|8.1| ) of this theory evaluates to zero. Performing the 
necessary subtraction we find 

r^black brane-sub ^ | ^_ ((2^T)'^+1 - l) 5^, + (d + 1) {{2'kTY+H ^t^.t - d^e^.e)] ■ (E.4) 

In particular, the energy density e, free energy density / and pressure P of the black brane 
(the deconfined phase) are given as functions of temperature by 

e=|(d(27rT)'^+i + l); P = -/ = | ((27rT)^+i - l) . (E.5) 

Appendix F. Numerical Construction of the Domain Wall Solution 
F.l. Equations and Constraints 

Working in units with LF'a' = 1, the metric we require is Einstein and must solve 
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(F.l) 



In our coordinate system (8^) the three equations from the rr, 66 and VaVa components 
manifest the eUiptic nature of the static Einstein equations, where the second derivatives 
are simply the flat two dimensional {x, 2/)-Laplacian, giving coupled second order equations 
for A, B,C. The equation for A is 

V^A ={d+ \)A^^ - [{d - 2) VC + ■ VA, (F.2) 

and the equation for B is identical after the substitution A B above, since in the static 
context there is no physical difference between the Euclidean time circle and the compact 
spatial circle 6. The equation for C is of a similar form, 

V^C = + l)e2" - f^d - 2) VC + ^ + ^) ■ VC. (F.3) 

The last elliptic (i.e. Laplacian second derivatives) equation can be found for the final 
function D from the xx + yy component of ([F.l| ) , 

(d+l)(rf-2) 2D , ,,f{d-3)^^ , VA , VB\ , 1 



V^D = q ^e- + id-2) [^VC + ^ + ^ J • VC + —VA • VB. 

(F.4) 

The remaining Einstein equations, namely the xy and the xx — yy components of ( [F.l| ), do 
not have elliptic second derivatives (but rather dxdy and — dy respectively) and should 
best be thought of as constraints associated with the gauge-fixing performed in (|8.6| ). Let 
us write them as 

a = Vdetg^ {R'y) = 0, 

1 , (F.5) 

The two non-trivial contracted Bianchi identities then relate a and /3 as a very elegant 
Cauchy-Riemann problem, 

d^a + dy(3 = 0, 

(F.6) 

dya - d^P = 0, 

provided that the above elliptic equations, (|F.2| ), ( [F.3|) and ( [F.4|) , are satisfied. 

Our strategy in finding solutions will be to solve the elliptic equations, subject to 
boundary data that ensures that the Cauchy-Riemann problem ( [l''.6| ) only has the trivial 
solution a = (3 = 0, and hence the constraint equations ( [t'^.S] ) will be satisfied. The natural 
data to take for these elliptic equations is to specify a Dirichlet or Neumann (or mixed) 
condition on each boundary for each metric function. 
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F.2. The IR Boundary Conditions 

We now consider the boundary conditions on the IR coordinate boundary where the 
r or ^ circle shrinks. Since these two circles are identical from the point of view of the 
equations, let us choose to discuss the shrinking r circle at y = (the other case being 
identical with x ^ y and A ^ B). 

There are two sources of data for the functions A, B, C, D. Firstly, the shrinking cycles 
imply regular singular behavior in the elliptic equations, and our choice that these metric 
components be finite singles out certain boundary conditions. Secondly, the constraint 
equations a = 0, (3 = also lead to boundary conditions. 

Firstly, consider the boundary conditions from the elliptic equations at y = 0. By 
construction we impose the Dirichlet data A = 0. When A 0, some terms in the 
'sources' for these elliptic equations ([F.2|) , (|F.3|) and (|F.4| ) are potentially singular, as they 
go as 1/A. Regularity therefore imposes Neumann boundary conditions on B and C, 

dyB = 0, 



(F.7) 

ai y = 0. Secondly, consider the constraints at y = 0. These give, 

(F.8) 



dyC = 0, 



a : (e-^'dyA) = 0, 



(3 : dyD = 0. 

We may immediately recognize the a = constraint as the zeroth law of black hole 
mechanics. Similarly, for the spatial circle shrinking on a; = 0, it ensures that the geometry 
will close smoothly for one value of the circle coordinate period. 

Note that these boundary conditions are derived from the equations of motion, and 
are totally independent of whether we work in the Euclidean or Lorentzian setting. 

Solving the elliptic equations in an elliptic manner we require one piece of data on 
each boundary for each metric function. Thus we see that we now have an over-complete 
set of boundary conditions for A, B,C,D in the IR. We therefore choose to impose only 
the P = constraint, and hence obtain the rather simple Neumann and Dirichlet data. 



(F.9) 



A 


= 0, 


dyB 


= 0, 


dyC 


= 0, 


dyD 


= 0, 
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at 2/ = 0, and similarly for the coordinate boundary at a; = the same with x ^ y and 
A ^ B. We leave the a = constraint to be imposed via the constraint system, by 
choosing data appropriately on the other boundaries. 

An important point is that at the phase interface x = y = both constraints are 
satisfied automatically due to the factor of AB occurring in ^/det g^^. The two constraints 
take the form, 

a = A ax=0 + B ay=o + AB aremainder, 

(F.IO) 

13 — A I3x=0 + B Py=0 + AB ^remainder, 

where ai, (3i are linear functions in derivatives of A and B, but are not explicit functions 
of A and B themselves. The terms ay=0: /3y=o give the constraints ( [l'\8| ) on the y = 
shrinking time circle boundary, but do not contribute to the shrinking space circle boundary 
a: = due to the prefactor of B (omitted in (|F.8| )). Likewise aa;=07 /3x=o give the constraints 
for the shrinking circle boundary x = 0, and don't contribute to the y = boundary. The 
final terms cxremainder, (^remainder Contribute to neither boundary as AB is zero on both, 
and also at the interface. Hence no additional data is required to satisfy the constraint 
system at the interface, beyond ( [F.8|) (and its symmetric counterpart at x = 0). 



F. 3. The UV Boundary Conditions 

The UV boundary resides at the locus / = where we demand that the metric looks 
like the boundary of AdS space, with A, S, e*^, ~ 1// (assuming that / has a first order 
zero). In order to impose this boundary condition, near the UV boundary we write ( |8.6| ) 
in the form 

ds^ = ^ {T^dr^ + U^d6^ + e^^drl + e^^lV/l^ (rfx^ + dy^)) (F.ll) 

and we take the new metric variables T, f/, S, R (which are trivially related to A, S, C, D 
via /) to be finite there. Coupled with the form above, we also make the further choice 
that the function / obeys the Laplace equation 

(a2 + aJ)/ = o, (F.12) 

and is defined in a finite neighborhood of the UV boundary. Note that this is consistent 
with the fact that for large x we expect / oc (y — c), and / oc (x — c) for large y. The 
reason for choosing / to be harmonic is that then we can easily map ( [F.ll[) to 



ds^ = -1 (T^dT^ + U^d6^ + e^^drl + e^^ {df + dg^)) (F.13) 
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where the metric components are now written as functions of f{x,y) and g{x,y), with 
the coordinate g defined by the Cauchy-Riemann relations, dxf = dyg and dyf = —d^g, 
so that the new coordinates /, g are conformaUy related to x, y. Since / = at the UV 
boundary, it can be thought of as the asymptotic bulk radial coordinate (like u earlier), 
and g is naturally then the coordinate along the boundary at this zero of /. 

In the coordinates ( [F.13|) it is easy to see that AdS space is simply given by 
T, f/, e"^ =constant and = 1, since then the metric (JFTT^) is just AdS written in the 



usual Poincare coordinates. Hence, T, t/, going to constants and ^ 1 near / = 
will reproduce asymptotically AdS space. Reversing the argument, starting with AdS in 
Poincare coordinates near the boundary and choosing our coordinate system (|8.6| ) and the 
position of the UV boundary to be fixed on some curve /(x, y) = 0, then the most general 
form that the asymptotically AdS region can take is with T, [/, =constant and 

= 1, for some harmonic function / with the correct zero locus. 
We now make the further choice that 

T,U,e^,e^^l (F.14) 

as / — > 0, for some given /, so that the metric of the dual field theory, in the conformal 
frame defined by /, is just 

dsl^^ = dr^ + de^ + drl + dg^. (F.15) 

This choice is always possible, as we may independently rescale the r, 9 and coordinates 
without effecting the Einstein equations. It is this choice (together with the choice of UV 
coordinate boundary location) that will set the actual values of /3, P required to smoothly 
close the manifold in the IR. Note that we have chosen the asymptotic values of T, U to 
be the same in order to preserve the natural ^2 symmetry of the solution. 

Suppose that for some metric A, B, C, D we choose a particular / satisfying (|F.12| ) 
such that the behavior (|F.14|) is true. Now let us consider a different /, say /', which 



shares the same locus for its zero, and also satisfies (IF. 121). Then the induced T' U' S' R' 



using /' will not satisfy ( |F.14 ), but rather will go as 



e^' ^ 1. 



|V/| ' (F.16) 



48 



This gives a boundary metric in the conformal frame defined by /' which is 



ds'i^^ = (j){hf {dr^ + de^ + drl + dh^) (F.17) 

with dg' = (f) dh. Hence, different choices of /, together with the condition ( |t'\14|) , pick 
out conformaUy related boundary metrics, where the normal gradient of / to its zero 
determines the conformal scaling. 

Thus we could pick a particular / and implement the freedom of the conformal factor 
in the asymptotic behavior directly (as in ([F.16| )). However this considerably complicates 
the asymptotic behavior of the metric functions. Instead, we do not pre-determine /, but 
impose the choice ([F.14| ) and then this scalar degree of freedom is just the normal gradient 
of / to its zero. 

For any / satisfying ( [F.12|) , the 'elliptic' Einstein equations for T, t/, e*^, given by 
( [F.2|) , (|F.3| ) and ([F.4|) , then yield the following simple behavior near / = : 



T = 1 + f^+H^g) + /'^+V,+2(^7) + /'+V,+3(^7) + ■ ■ ■ , 

U = l + f+'^ig) + f+'u.+^ig) + f+'u.+sig) + ■ ■ ■ , 

1 (F.18) 

{d-2) 

e« = 1 + + + f^^ra+M + • " " , 

where again g is b. coordinate normal to / which measures the position along the locus of 
the zero of /. The four leading ones in these expansions should be likened to Dirichlet 
data, and as discussed above, guarantee that this geometry is asymptotic to AdS. The 
four undetermined finite functions 5, 7,p, a should be regarded as Neumann data for the 
elliptic equations. Since we fix the Dirichlet data, we are not at liberty to also fix these 
functions when solving the system in an elliptic manner, and these will be determined 



Choosing the three leading (Dirichlet) terms in the asymptotic expansions of r, U, to be 
finite functions of g (rather than one as in ( 1F.18D ) represents a 'non-normahzable' deformation of 
the asymptotic region of the space away from AdS, and will give a non-flat boundary metric. The 
fourth piece of Dirichlet data is implicit in the finiteness of R at the UV boundary. Whilst ^ 1 
if R is finite, there is an additional possible non-normalizable behavior, seen most easily in the 
linear theory where R r{g) log(/) is also possible. This can simply be understood as arising 
from a conformal coordinate transformation moving the coordinate position of the UV boundary, 
and hence being singular in R. This function r{g) which we have implicitly set to zero in (F.18) 
is the fourth and last piece of Dirichlet data for the functions T, [/, S, R at the UV boundary. 
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implicitly by the IR boundary conditions. The functions {r, m, r, s}i for i > d + 2 are terms 
in the power-series expansion of the metric in / which are determined in terms of these 
leading 'Dirichlet' and 'Neumann' coefficients as in the usual Fefferman-Graham expansion 



Using these expansions we may compute the boundary stress tensor as in appendix 
E, with components given by 

T^r = {d+mg) + p{g), 
Tee = id+lh{g) + p{g), 

Tr.r. = - (id + img) + -fig)) + {d + 4)p{g)) , ^^'^^^ 

Tgg = id + 2)pig) 

in an appropriate normalization. Note that, as required, this is traceless. Recall that we 
should subtract the confining vacuum stress tensor in order to obtain the standard dual 
field theory stress tensor. In the fiat boundary metric, the conservation equation of the 
stress-energy tensor gives dgTgg = 0, and this is crucial in order to understand the pressure 
balance argument which singles out a particular temperature thermal bubble interior. In 
the gravity context this arises as a result of the equation of motion of G-^g, as we now 
discuss. 

So far we have only discussed the behavior of the 'elliptic' equations. Now we evaluate 
the constraint equations 5 = 0, (3 = 0, in the new coordinates /, g : 



a = ^/detg^ [R^g) = 0, 

~ 1 / . \ (F-20) 

(3 = -y^detg^^ [^Rff - R%j = 0. 

We note that 5, (3 are simply linearly related to a, (3, the coefficients depending on the 
(non-singular) coordinate transformation between the /, g and x, y coordinates. Using the 
expansion (|F.18| ) we find that both these constraints are finite near / = 0, although neither 
is automatically zero, with 

---'-^^^M^oif), 

(3={d + 2)a{g)+0{f). 
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We focus on the 5 = constraint, which is particularly elegant as it implies that p{g) is a 
constant, say k, and hence contours of / and R are aligned asymptotically. Thus it implies 
that near the zero of /, 

R{f,9) = kf''+' + 0{f+'). (F.22) 

Furthermore, this obviously guarantees that Tgg =constant in the boundary stress tensor 
(and will be zero with a suitable vacuum subtraction). Note that the actual value of k is 
unphysical since in the definition of / via ( [F.ll|) we see the physical metric is invariant un- 
der rescalings of /, while appropriately rescaling the t,9 and coordinates too. Therefore 
we may choose any value for k that is convenient. 

F.4- Solving the Constraint System 

At large x, y we expect our solution to return to the homogeneous black brane ( |875| ) 
or to the AdS soliton (|8.4| ) (this should occur exponentially fast in proper distance, since 
the effective dual theory has a mass gap). Since the a = constraint is satisfied for 
a homogeneous solution (depending only on x or y), at large x,y where the solution 
should become homogeneous it will quickly go to zero. By the Cauchy-Riemann constraint 
relations ( [F.6|) /3 must then go to a constant, and since we are actually imposing /3 = on 
the IR coordinate boundary, it will also vanish. Hence, the constraints will be satisfied in 
these asymptotic regions. In the IR we have explicitly imposed jS = on the axes a; = 
and y = 0. This is still not sufficient data to solve the constraint problem. 

In order to solve the system we could try to impose a in the IR, but as we discussed 
above, the elliptic system is already completely determined there. Instead, we may impose 
the vanishing of some linear combination of a and /3 in the UV. With /3 = in the IR, and a 
linear combination being zero in the UV, combined with the correct large x, y asymptotics, 
this is sufficient to guarantee that the solution of the Cauchy-Riemann system ( [b'.GD is such 
that both a and /3 vanish everywhere. 

From the discussion above, the natural combination of a and (3 to impose is exactly 
5, coming from the fg component of the Einstein equations in the UV, which implies the 
relation ( [l'\22| ) . Now we already have Dirichlet data for T, t/, S, R in the UV, determined 
by leading behavior in (|F.18| ), and we cannot also fix the Neumann terms 5, 7,p, cr. The 



only undetermined quantity is the choice of our harmonic function /. Thus, we impose 
the constraint 5 = in the UV by choosing / to have the correct normal gradient solving 
( [b'.22D on its zero locus, so 

|V(/'^+^)| = ]-\VR\ (F.23) 



k' 
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when / = 0. Note that since the Laplace equation is eUiptic, this data does not determine 
/ globally in our coordinate domain, but it is sufficient to determine it exponentially well 
near the UV boundary, and we only require / in the immediate vicinity of this boundary 
(i.e. we really only require its normal derivative) to define the asymptotic behavior of 
A, S, C, D. The behavior (|F.18| ) will now ensure that the constraint 5 = is satisfied. 
As mentioned above, the constant of proportionality, k, is unphysical and can be chosen 
arbitrarily. 



F.5. Summary 

Our complete system, in either Euclidean or Lorentzian signature, is then the elliptic 
equations ([F.2| ),( p\3D ,( [F.4| ) and ([F.12| ), together with their boundary conditions ( [F.9|) , 
( [b'.14[ ) and ( [1''.2^ ). We see explicitly that, as indicated in section 8, the 'Z2 symmetry ( P78| ) 
is manifest in the system. Furthermore, for fixed boundary locations, and once we have 
fixed the arbitrary scaling symmetry of the r, 9 and coordinates (through the combined 
choices of c, k and the choice ( [F.14|) ), there is no further data to be specified. Whilst 
we have not rigourously demonstrated that this system is elliptic, this analysis strongly 
suggests that the solution is unique, having no additional parameters. 



Appendix G. Numerical Details 

We define our coordinate domain using an auxiliary function h. We set /i = 1 on the 
lines a; = and y = 0, and h = —0.6 at x > 1, y = 1 and y > 1, x = 1. We solve a Laplace 
equation for h in the interior of this region. Using this solution we then determine the 
coordinate position of the UV boundary at /i = where we will demand / = 0, the IR 
boundaries at a; = and y = 0, and hence the coordinate domain for our problem. Note 
this gives the required ^2 symmetry of the UV boundary location, and for these choices, 
at large x the UV boundary tends to y = 5/8. 

Furthermore we also determine an interior 'boundary' at h{x, y) = ho, which satisfies 
< /io"*"^ << 1 for computational purposes. To the IR of this interior boundary {h > ho) 
we use the metric variables A, S, C, D when solving the elliptic equations, and to the UV 
{h < ho) we use T, C7, S, R which allows us to encode the correct UV behavior, without 
requiring large gradients that reduce accuracy and stability. The function / is now only 
defined in this UV region < h < ho, where we must be able to translate between 



52 



A, B, C, D and T, U, R. The initial guess for / is then taken to be h, and sensible initial 
guesses for the metric functions A, B, C, D and T, U, S, R are then made (using h). 

Due to the symmetry (|8.8| ) , we only implement the region y < x and use the symmetry 
to impose boundary conditions at x = y. In practice we truncate the asymptotic region 
at large x = Ltmnc and simply require that normal gradients of all functions vanish 
there (i.e. Neumann boundary conditions). We also truncate the UV boundary to be at 
/ = e, rather than / = 0, since the behavior of various terms in the elliptic equations 
are potentially rather singular there. We check that when the position of the large x, y 
boundaries is large compared to the bulk curvature length this truncation does not influence 
the solution. Similarly, when e is small, the difference between the solution for finite e and 
the extrapolated e = solution is also negligible. 

The elliptic equations are second order finite differenced. The metric functions are 
relaxed using Gauss-Seidel iteration with the boundary conditions discussed earlier. Some 
under-relaxation is required depending on the resolution used. For some number of itera- 
tions / is not updated from its initial guess, to allow the metric components to settle to 
have sensible asymptotic behavior. Once this has occurred, / is also relaxed in the small 
UV region by the same method. 

The value of / at the interior boundary h = ho is required to solve R = kf^^^ there 
for some value of k. We solve this condition by updating / on this interior boundary by 
5f = uj{R/kf'^^^ — 1), where 5f is the update for / after each Gauss-Seidel step and ui 
is a small number (say ~ 10~^ for the resolution used to produce the data shown here). 
This very slow update is required as the UV boundary is rather unstable numerically due 
to the singular nature of the elliptic equations there. 

Ideally we only require the normal gradient behavior of / and thus really require only 
a tiny h^, so the UV region where / is defined is very small. In practice, the singular 
behavior of the equations near the UV boundary mean /iq cannot be taken too small at 
a given resolution. However, we checked that for the values and resolutions used, the 
solutions were insensitive to the choice of /iq, provided it is small enough. 

Results shown in the paper use the resolution 320 x 80 in x, y (recall we only use the 
region y < x) with Ltmnc = 4, e = 0.05 and /iq = 0.6 for both (i = 3, 4. The values of k 
were chosen for convenience to be 0.35 and 0.3 for d = 3,4 respectively. 
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Figure 12: Contour plots of A, S, e*-" , in tlie x, ?/ plane for d = A. 
We remind the reader that the explicit program code used to generate the solutions 



may be downloaded at [http : / / schwinger . harvard . edu/^wiseman/IRblackholes/ 



The procedure works very well, converging to a unique solution independent of initial 
guesses. We now give the results in more detail. The earlier figures 5 and 6 showed surface 
plots of the functions / and A, S, e*^, for d = A. In figure 12 we show contour plots of 
the functions A, S, e*-" , for d = A. The functions appear qualitatively similar for (i = 3, 
and can not be distinguished from these d = A contour plots by eye. 

Let us now check that the 5 = constraint is correctly imposed in the UV. In figure 
13 we plot, for (i = 4, contours of R and / near the UV boundary, and indeed they agree 
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Figure 13: Plot showing contours of constant / (red, dashed) and R (black, 
solid) for d = 4 near the UV boundary where / is defined. These are aligned, 
indicating that the constraint S = is satisfied near the UV boundary. 
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Figure 14: Plot showing violation of the a = constraint in the IR. The 
proper gradient of the vanishing circle size, k, is plotted against position on 
the boundary x for d = 3 (left) and d = 4 (right), normalized to its value as 
X — > 00. The a = constraint implies that this should be independent of x. 
Recall that we impose the constraints /3 = in the IR and 5 = in the UV. 
Here we see that a = is true to high accuracy in the IR, implying we have 
indeed solved the Cauchy-Riemann constraint system. 

very well. One obtains a very similar plot for d = 3. 

We directly impose the /3 = constraint in the IR and the 5 = constraint in the 
UV. Thus we do not impose the a = constraint in the IR and this gives a rather physical 
test of how well we solve the constraint system ( [F.6| ) by plotting the proper gradient of 



A at the IR boundary y = 0, i.e. k{x) = e dyA\y^o. We recall from (II^'.SI ) that the 



a = constraint in the IR implies that k(x) should in fact be a constant, and 2n / k gives 
the coordinate period of the two circles /3, (3. In figure 14 we plot k{x) for the d = 3 and 
d = 4 solutions, and we see that they are both beautifully constant indeed, refiecting the 
fact that we have solved the Cauchy-Riemann constraint system to a high accuracy for the 
resolution used. 
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